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Abstract 

A numerical algorithm for explicitly computing the spectrum of the 
Laplace-Beltrami operator on Calabi-Yau threefolds is presented. The 
requisite Ricci-flat metrics are calculated using a method introduced in 
previous papers. To illustrate our algorithm, the eigenvalues and eigen- 
functions of the Laplacian are computed numerically on two different quin- 
tic hypersurfaces, some Z5 x Z5 quotients of quintics, and the Calabi-Yau 
threefold with Z3 x Z3 fundamental group of a heterotic standard model. 
The multiplicities of the eigenvalues are explained in detail in terms of the 
irreducible representations of the finite isometry groups of the threefolds. 
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1 Introduction 



A central problem of string theory is to find compactifications whose low-energy ef- 
fective action reproduces the standard model of elementary particle physics. One 
of the most promising candidates for this task is the compactification of heterotic 
string theory on a Calabi-Yau manifold [1]. In particular, the so-called "non-standard 
embedding" of Eg x Eg heterotic strings has been a very fruitful approach to string 
phenomenology [21 S 01 El El E E] • 

For a number of reasons, the most successful models of this type to date are 
based on non-simply connected Calabi-Yau threefolds. These manifolds admit discrete 
Wilson lines which, together with a non-flat vector bundle, play an important role in 
breaking the heterotic Eg gauge theory down to the standard model [9], UHl HH [121 [HI 
[HI, [151 [16] . Ill addition, they project out many unwanted fields which would otherwise 
give rise to exotic matter representations and/or additional replicas of standard model 
fields. In particular, one can use this mechanism to solve the doublet-triplet splitting 
problem [TTIIIH]- Finally, the non-simply connected threefolds have many fewer moduli 
as compared to their simply connected covering spaces [19] . In recent work [201 [211 
[22l [23] . three generation models with a variety of desirable features were introduced. 
These are based on a certain quotients of a Schoen Calabi-Yau threefold, yielding a 
non-simply connected Calabi-Yau manifold. 

The ultimate goal is to compute all of the observable quantities of particle physics, 
in particular gauge and Yukawa couplings, from the microscopic physics of string 
theory [211 [2Sl [26l [2Z]- There are many issues which must be addressed to achieve 
this goal. Physical Yukawa couplings, for example, depend on both coefficients in the 
superpotential and the explicit form of the Kahler potential. In a very limited number 
of specific geometries [211 [23 [23 [30] , the former can be computed using sophisticated 
methods of algebraic geometry, topological string theory and the like. For the latter, 
one is usually limited to the qualitative statement that a coefficient is "expected to be 
of order one" . Improving our computational abilities and extending these calculations 
to non-standard embedding has been an outstanding problem [1]. 

Recently [SH [32], a plan has been outlined to analyze these problems numerically, 
at least in the classical limit. The essential point is that, today, there are good enough 
algorithms and fast enough computers to calculate Ricci-flat metrics and to solve the 
hermitian Yang-Mills equation for the gauge connection directly. Given this data, one 
can then find the correctly normalized zero modes of fields, determine the coefficients 
in the superpotential and compute the explicit form of the Kahler potential. Some 
progress in this direction was made in [3ll [321 [331 [311 [3H] and also [33 [33, [3H] • Making 
effective use of symmetries [5^ BU] , one can significantly improve the computational 
procedure to find Calabi-Yau metrics and further extend it to non-simply connected 
manifolds. In this work, we take one step further in the numerical approach to string 
theory compactification and present an explicit algorithm to numerically solve for the 
eigenvalues and eigenfunctions of the scalar Laplace operator. We use as one of the 
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inputs the Calabi-Yau metrics computed using the techniques developed in [ID]. 

We start, in [Section 31 by discussing the general idea of the method and list the 
key steps of our algorithm. This algorithm is then applied to the simplest compact 
threefold, the projective space Pi This threefold is, of course, not a Calabi-Yau man- 
ifold. However, in has the advantage of being one of the few manifolds where the 
Laplace equation can be solved analytically. We compare the numerical results of 
this computation with the analytical solution in order to verify that our implementa- 
tion is correct and to understand the sources of numerical errors. We note that the 
multiplicities of the approximate eigenvalues are determined by the dimensions of cor- 
responding irreducible representations of the symmetry group of the projective space, 
as expected from the analytical solution. We conclude the section by investigating the 
asymptotic behavior of the numerical solution and comparing it with Weyl's formula. 

Having gone through this illustrative example, we apply our numerical procedure 
to Calabi-Yau quintic threefolds in [Section 4[ The eigenvalues and eigenfunctions are 
explicitly computed for both a quintic at a random point in moduli space as well as 
for the Fermat quintic. We can again explain the multiplicities of eigenvalues on the 
Fermat quintic as arising from its enhanced symmetry; here, however, being a finite 
isometry group. The asymptotics of the numerical solution is verified using Weyl's 
formula. Note that the eigenvalues and eigenfunctions are not known analytically 
in the Calabi-Yau case, so our numerical algorithm is essential for their calculation. 
Recently, Donaldson has proposed a different algorithm to solve for the spectrunfl of 
the scalar Laplacian. At the end of the section, we use it to numerically compute 
the eigenvalues and eigenfunctions on a random quintic and on the Fermat quintic 
and compare these to our results. In [Section 5| we consider non-simply connected 
Calabi-Yau manifolds, namely Z5 x Z5 quotients of certain quintic threefolds. The 
eigenvalues and eigenfunctions of the Laplacian are numerically computed using our 
algorithm, exploiting the Hironaka decomposition discussed in our previous paper ^0]. 
In this case, the multiplicities of the eigenvalues are determined by finite "pseudo- 
symmetries" [H]. We work out the necessary representation theory and again find 
perfect agreement with the multiplicities predicted by our numerical computation of 
the eigenvalues. We conclude this section by studying the moduli dependence of the 
eigenvalues for a one-parameter families of quintic quotients. 

In [Section 6l we apply this machinery to the case of a certain Z3 x Z3 quotient of 
a Schoen threefold [l2l US] . This is the Calabi-Yau threefold underlying the heterotic 
standard model constructed in [211 ISH [23] . The essential new feature is the existence 
of non-trivial Kahler moduli, not just the overall volume of the threefold as in all 
previous sections. As an explicit example, we numerically compute the eigenvalues 
of the Laplacian at two different points in the Kahler moduli space, corresponding 
to distinct "angular" directions in the Kahler cone. The group representation theory 
associated with the covering space and the quotient is discussed. 

^The spectrum of an operator is the set of eigenvalues. 
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We conclude in lSection 71 by considering some physical applications of the eigenval- 
ues of the scalar Laplacian on a Calabi-Yau threefold. In particular, we consider string 
compactifications on these backgrounds and study the effect of the massive Kaluza- 
Klein modes on the static gravitational potential in four-dimensions. We compute 
this potential in the case of the Fermat quintic, and explicitly show how the potential 
changes as the radial distance approaches, and passes through, the compactification 
scale. We then give a geometrical interpretation to the eigenvalue of the first excited 
state in terms of the diameter of the Calabi-Yau manifold. Inverting this relationship 
allows us to calculate the "shape" of the Calabi-Yau threefold from the numerical 
knowledge of its first non-trivial eigenvalue. 

Additional information is provided in three appendices. We explicitly determine 
the first massive eigenvalue for the Laplacian P"^ in Appendix A Some technical 



aspects of semidirect products, which are useful in understanding [Section 4l are dis- 
cussed in Appendix B Finally, in Appendix C we explain a modification of Donald- 
son's algorithm for the numerical computation of Calabi-Yau metrics on quotients, 
which is used [Section 5l 



2 Solving the Laplace Equation 

Consider any d-dimensional, real manifold X. We will only be interested in closed 
manifolds; that is, compact and without boundary. Given a Riemannian metricH g^^ 
on X, the Laplace-Beltrami operator A is defined as 

A = -^d^ig^-'^^d,) = -5d= - * d* d, (1) 

where g = detg^y. Since this acts on functions, A is also called the scalar Laplace 
operator. We will always consider the functions to be complex-valued. Since A 
commutes with complex conjugation, the scalar Laplacian acting on real functions 
would essentially be the same. 

An important question is to determine the corresponding eigenvalues A and the 
eigenf unctions (j) defined by 

A0 = A(/>. (2) 

As is well-known, the Laplace operator is hermitian. Due to the last equality in 
eq. ([T]), all eigenvalues are real and non-negative. The goal of this paper is to find the 
eigenvalues and eigenf unctions of the scalar Laplace operator on specific manifolds X 
with metrics g^^. 

Since X is compact, the eigenvalues of the Laplace operator will be discrete. Let 
us specify the n-th eigenvalue by A„. Symmetries of the underlying manifold will, in 
general, cause A^ to be degenerate; that is, to have multiple eigenfunctions. We denote 

^We denote the real coordinate indices by n, v, .... 
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by /i„ the multiplicity at level n. Each eigenvalue depends on the total volume of the 
manifold. To see this, consider a linear rescaling of distances; that is, let g^i, h-> p^g^^. 
Clearly, 

Vol (p'^g^.u) = p'^ Vol [g^^) , A„ {p'^g^^y) = p~'^Xn {gf^u) ■ (3) 
Therefore, each eigenvalue scales as 

A„ ~ Vol"^ . (4) 

In the following, we will always normalize the volume to unity when computing eigen- 
values. 

Now consider the linear space of complex-valued functions on X and define an 
inner product by 

(e|/)= / ef^d\ e,/eC°°(X,C). (5) 
Jx 

Let {fa} he an arbitrary basis of the space of complex functions. For reasons to become 
clear later on, we will primarily be working with bases that are not orthonormal with 
respect to the inner product eq. ([5]). Be that as it may, for any complex function e 
one can always find a function e so that 

e = $^/a(/a|e). (6) 

a 

Given the basis of functions {fa}, the matrix elements A^;, of the Laplace operator 



are 



(/„|A|/,>= [ faAf, ^d'x = - [ fad*dft= [ {dfa\df,) 

Jx Jx Jx 



Thus far, we have considered arbitrary ci- dimensional, real manifolds X and any 
Riemannian metric gfj,iy. Henceforth, however, we restrict our attention to even dimen- 
sional manifolds that admit a complex structure preserved by the metric. That is, 
we will assume that X is a D = ^-dimensional complex manifold with an hermitiai]j 
metricQ gij defined by 

g^^ dx^ (g) dx" = ^gi-j ( dz' ® dz^ + dz^ ® dz') . (8) 
With X so restricted, it follows that 

g'^'^d^fa djh = 2g'^ (d-Ja d,f, + dja d-Jh) (9) 

■^In particular, Kahler metrics are hermitian. 

^We denote the holomorphic and anti-holomorphic indices hy i, t, j, j, .... 
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and, hence, 

Aab = 2 g'^ (d-Ja d,h + dja O-Jb) det{g) {^fH dz"- A d/. (10) 

r=l 

Using this and eq. (E]) for each eigenfunction (j)n,i, eq. ([2]) becomes 

J2 {fa\A\h){Mk,) = Yl Ufa\fb){fb\4>n,), ^ = 1, • • • , (H) 

b b 

Thus, in the basis {fa}, solving the Laplace eigenvalue equation is equivalent to the 
generalized eigenvalue problem for the infinite dimensional matrix Aab, where the 
matrix {fa\fb) indicates the "non-orthogonality" of our basis with respect to inner 
product eq. ([5]). 

In general, very little known about the exact eigenvalues and eigenf unctions of 
the scalar Laplace operator on a closed Riemannian manifold X, including those 
that are complex manifolds with hermitian metrics. The universal exception are the 
zero modes, where the multiplicity has a cohomological interpretation. Specifically, 
the solutions to = are precisely the locally constant functions and, hence, the 
multiplicity of the zero eigenvalue is 

fio{X) = h%XX) = \MX)\, (12) 

the number of connected components of X. Furthermore, on symmetric spaces G/H 
one can completely determine the spectrum of the Laplace operator in terms of the 
representation theory of the Lie groups G and H. Indeed, in the next section we 
will discuss one such example in detail. However, in general, and certainly for proper 
Calabi-Yau threefolds, exact solutions of = X(f) are unknown and one must employ 
numerical methods to determine the eigenvalues and eigenf unctions. The purpose 
of this paper is to present such a numerical method, and to use it to determine 
the spectrum of A on physically relevant complex manifolds. Loosely speaking, the 
algorithm is as follows. 

First, we specify the complex manifold X of interest as well as an explicit hermitian 
metric. For Kahler manifolds, the Fubini-Study metric can always be constructed. 
However, this metric is never Ricci-flat. To calculate the Ricci-flat Calabi-Yau metric, 
one can use the algorithm presented in [311 ES] and extended in [ID] . This allows a 
numerical computation of the Calabi-Yau metric to any desired accuracy. Giving the 
explicit metric completely determines the Laplace operator A. Having done that, we 
specify a countably infinite set {fa} that spans the space of complex functions. One 
can now calculate any matrix element Aab = {fa\A\fb) and coefficient {fa\fb) using 
the scalar product specified in eq. ([5]) and evaluated using numerical integration over 
X. As mentioned above, the most convenient basis of functions {fa} will not be 
orthonormal. Clearly, calculating the infinite dimensional matrices Aab and {fa\fb), 
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let alone solving for the infinite number of eigenvalues and eigenfunctions, is not 
possible. Instead, we greatly simplify the problem by choosing a finite subset of 
slowly-varying functions as an approximate basis. For simplicity of notation, let us 
take {fa\a = 1, . . . , A;} to be our approximating basis. The k x k matrices {Aab)i<a,b<k 
and {fa\fb)i<a,b<k are then finite dimensional and one can numerically solve eq. f|TT]) 
for the approximate eigenvalues and eigenfunctions. It is important to note that 
this procedure generically violates any underlying symmetries of the manifold and, 
hence, each eigenvalue will be non-degenerate. Finally, we successively improve the 
accuracy of the approximation in two ways: 1) for fixed k the numerical integration of 
the matrix elements is improved by summing over more points and 2) we increase the 
dimension k of the truncated space of functions. In the limit where both the numerical 
integration becomes exact and where A; — > oo, the approximate eigenvalues A„ and 
eigenfunctions 0„ converge to the exact eigenvalues Xm and eigenfunctions (j)rn,i with 
multiplicity Inspired by our work on Calabi-Yau threefolds, this algorithm to 
compute the spectrum of the Laplacian was recently applied to elliptic curves in [H]. 

3 The Spectrum of A on 

In this section, we use our numerical method to compute the eigenvalues and eigen- 
functions of A on the complex projective threefold 

p3 = syu{l) = SU{A) /S{U{3) X f/(l)) (13) 

with a Kahler metric proportional to the Fubini-Study metric, rescaled so that the 
total volume is unity. As mentioned above, since this is a symmetric space of the form 
G/H, the equation A0 = A0 can be solved analytically. The results were presented 
in ^5\. Therefore, although is not a phenomenologically realistic string vacuum, 
it is an instructive first example since we can check our numerical algorithm against 
the exact eigenvalues and eigenfunctions. Note that, in this case, the metric is known 
analytically and does not need to be determined numerically. 

3.1 Analytic Results 

Let us begin by reviewing the known analytic results [15]. First, recall the Fubini- 
Study metric is given by gfj^ = didjKps with 

Kfs{z, ^) = ^ In (kor + kiT + 1^21' + Nsr) • (14) 
With respect to this metric the volume of is 

Volps(P=^)= / det{g.^d'x= / = ^, (15) 
Jp3 Jp3 3! b 
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where wps is the associated Kahler (1, l)-form. However, as discussed above, we find 
it convenient to choose the metric so as to give unit volume. It follows from eq. (fT^ 
and (1151) that one must rescale the Kahler potential to be 

K{z,z) = ^Kfs{z,z) = (^01' + + ^2!' + ksl')- (16) 

Then 

VoW(P=^) = l, (17) 

as desired. 

The complete set of eigenvalues of A on P'^ were found to be [IS] 

Air 

A„ = -=m(m + 3), m = 0,1,2,..., (18) 
V6 

where we determine the numerical coefficient, corresponding to our volume normal- 



ization, in Appendix A Furthermore, it was shown in [15] that the multiplicity of the 



m-th eigenvalue is 



lim=\ - J =-(m + l)^(m + 2)^(2m + 3). (19) 

This result for the multiplicity has a straightforward interpretation. As is evident from 
the description of P^ in eq. f|T3l) . one can define an SU{A) action on our projective 
space. Thus the eigenstates of the Laplace operator eq. ([2]) carry representations of 
SU (4). In general, any representation of S'f/(4) is characterized by a three dimensional 
weight lattice. In particular, for each irreducible representation there exists a highest 
weight 

w = miWi + m2W2 + rn^ws, (20) 

where Wi, W2, and W3 are the fundamental weights and mi,m2,m-^ G Z>o. Starting 
with the highest weight, one can generate all the states of the irreducible representa- 
tion. It turns out that multiplicity eq. f|T9|) is precisely the dimension of the irreducible 
representation of S'f/(4) generated by the highest weight m{wi + W3) = (m, 0,m). 
Hence, the eigenspace associated with the m-th eigenvalue carries the irreducible 
representation (m, 0,m) of SU{4) for each non-negative integer m. For convenience, 
we list the low-lying eigenvalues and their corresponding multiplicities in ITable II 

The eigenfunctions of A on P^ = S'^/U{1) are the f/(l)-invariant spherical har- 
monics on S'^. In terms of homogeneous coordinates [20 : Zi : Z2 '■ Z3] on P^, the 
eigenfunctions can be realized as finite linear combinations of functions of the form@ 



degree monomial ) ( degree ks monomial 



(21) 



^We label the degree of the monomials here by to distinguish it from the degree kh of polyno- 
mials in Donaldson's algorithm. 
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m 



1 

2 

3 

4 

5 

6 

7 



1 

15 

84 
300 
825 
1911 
3920 
7344 





n - 27.662 
^ - 69.155 

727r , 
1127r 
1607r 

TT 

21677 

'W 

2807r 



124.48 
193.64 
276.62 
373.44 
484.09 



Table 1: Eigenvalues of A on ¥ . Each eigenvalue is listed with its multiplicity. 

One can show this as follows. Let 4 and 4 be the fundamental representations of 
SU{4). Algebraically, one can show that 



Sym'''^ 4 (g) Sym'''^ 4 = fl^ (m, 0,m 



(22) 



where (m, 0, m) are the irreducible representations of SU{4) defined above. Now note 
that C[2|fc^, the complex hnear space of degree- /c^ homogeneous polynomials in zq, 
zi, Z2, -Zs, naturally carries the Sym'^'* 4 reducible representation of SU{A). Similarly, 
C[l]fc_^ carries the Sym'^'''4 representation. Defining 



, Zi, Z2, Z^ , ^i, Z2, ^Slkj, 



(23) 



to be the space of functions spanned by the degree monomials, then it follows from 
eq. fl22p that one must have the decomposition 



5" k 



span 



m,l) • • • ; Ym,fim J ) 



} 



(24) 



m=0 



where /x^ = dim(m, 0,m). Note the importance of the S'f/(4)-invariant denominator, 
which ensures that the whole fraction is of homogeneous degree zero, that is, a function 
on 

To illustrate this decomposition, first consider the trivial case where = 0. 
Noting that /io = 1, eq. (!24l) yields 



^0,1 



(25) 
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corresponding to the trivial representation 1^ of SU (4) and the lowest eigenvalue Aq = 
0. Now, let krf) = 1. In this case /io = 1 and /ii = 15. It follows from eq. (IMl) that 
there must exist a basis of composed of the eigenf unctions of A in the 1 and 15 
irreducible representations of SU{4) respectively. This is indeed the case. We find 
that one such basis choice is 



^0,1 





2 + 






1^2 


? + 




2 


y3 1 




2 



(26) 



corresponding to the lowest eigenvalue Aq = 0, and 



^1,1 

01,3 
01,5 
01,7 
01,9 
^1,11 



^0^1^ 
^0^2^ 

^1^2^ 
2:1^3^ 
2:2^3, 



j=0 



E 
E 
E 
E 



3 

j=0 



3 

i=0 
3 

i=o 
3 

i=o 



3 



<A'l,2 
01,4 
01,6 
01,8 
'^1,10 
!'l,12 



^^l^^O^ 
^2^0^ 
2:3^0^ 
2:2^1^ 
^3^1^ 
^3^2, 



3 

i=o 



E, 
E, 
E, 
E, 
E, 
Ej=o 



3 

i=o 
3 

i=o 
3 

i=o 
3 

i=o 



2;, I 



(27) 



-'1,13 



■^1,14 



■^1,15 



2^1 2^1 — 2^0 ^0 
Z2Z2 — ZqZq 











)/t.U 


Zj\ 



corresponding to the first non-trivial eigenvalue Ai 



Note that we recover the 

v6 

constant eigenfunction for = through the cancellation of the numerator in eq. ( l26l) . 
This pattern, where one recovers all the lower eigenmodes through the factorization of 



the numerator in each representation by an appropriate power of ^ 
for arbitrary k^- In other words, there is a sequence of inclusions 



3 

j=o \ 



continues 



Note that 



{1} = 9^0 C :?i C 9^2 C ■■■ C C°"(P^C). 



(2^ 



(29) 



which, together with eq. ( !22ll . explains the multiplicities given in eq. ( fT9l) . 

Although a basis of 3"^^ composed of eigenfunctions of A would be the most natu- 
ral, there is no need to go through the exercise of decomposing the space into SU (4)- 



10 



irreducible representations. For numerical calculations, it is simpler to use the equiv- 
alent basis 

^■fc^ = span {/a I a = 0, ... , dim 5"^^ - l} 

^degree fc<)i, monomial j ^degree monomial j / (^^^kjpj 1" 

for any finite value of kfp, even though these functions are generically not themselves 
eigenfunctions of A. In the limit where — >• oo, the basis eq. (l30l) spans the complete 
space of eigenfunctions. 



3.2 Numerical Results 

Following the algorithm presented at the end of the [Section 2| we now numerically 
solve the eigenvalue problem for the scalar Laplace operator A on P^. Unlike more 
phenomenologically interesting Calabi-Yau threefolds, where one must numerically 
compute the Kahler metric using Donaldson's method [3ll [33l EO] , on the Kahler 
potential is given by eq. (fT6|) and, hence, the metric and A are known explicitly. This 
eliminates the need for the first few steps of our algorithm, greatly simplifying the cal- 
culations in this section. Furthermore, the SU (4) action on the eigenfunctions allows 
us to identify a complete basis for the space of complex functions in terms of mono- 
mials of the form eq. ( 1211) . Since we know the exact eigenvalues and eigenfunctions 
on this is an excellent venue for checking the numerical accuracy of the remaining 
steps in our algorithm as well as the correctness of our implementation. 

Given the metric, A and the complete basis of functions, the next step in our algo- 
rithm is to specify an approximating basis for the linear space of complex functions. 
This is easily accomplished by restricting to 

S^k, = span {/, I a = 0, ... , ftf )' - l}, (31) 

see eq. ( I30l) . for any finite value of k^. Next, we need to specify the volume measure 
in the integrals required to evaluate the matrix elements {fa\A\fb) and {fa\fb)- Each 
matrix element requires one integral over P^, as in eq. ([7]). The volume form is 
completely determined by the metric to be 

dVoW = (32) 

where ui is the Kahler (1, l)-form given by the Kahler potential eq. (fT6l) . Although P^ 
is simple enough to employ more elaborate techniques of integration, we will use the 
same numerical integration algorithm as with Calabi-Yau threefolds later on. That 
is, we approximate the integral by summing over random points. 



(33) 



11 



where / is an arbitrary function on The integration measure dVol in eq. ( !33l) is 
determined by the distribution of points. In other words, the random distribution of 
points must be chosen carefully in order to approximate the integral with our desired 
volume form dVoli<:. However, this can easily be done: simply pick the points in 
an SU (4)-uniform distribution. The corresponding integral measure is (up to overall 
scale) the unique S'?7(4)-invariant volume form, the Fubini-Study volume form. The 
normalization is fixed by our convention that Voli^(P^) = 1. 

The process of numerically evaluating integrals by summing over a finite number 
of points has one straightforward consequence. As discussed above, in the analytic 
solution the m-th eigenvalue Am is degenerate with multiplicity fim given in eq. ( fT9i) . 
The reason for the degeneracy is that the m-th eigenspace carries the (m, 0, m) highest 
weight representation of SU{A). However, even though the points have an SU{4)- 
uniform distribution, the simple fact that they are finite explicitly breaks the S'f/(4) 
symmetry. The consequence of this is that the degeneracy of each eigenvalue is com- 
pletely broken. It follows that in the numerical calculation, instead of one eigenvalue 
Xm with multiplicity /i^, one will find /i^ non-degenerate eigenvalues A„. Only in the 
limit that ^ oo will these converge to a single degenerate eigenvalue as 

Ao = Ao, . . . , A^„_i 

Ai,...,Ai5 = A^o, . . . , A^jo+^ji_i 

A165 • • • 5 A99 = A^o+^^, . . . , A^Q+^^+^2_i 



0, 

167r 

\ _ 407r 

^2 - 17g> 



Ao 
Ai 



(34) 



We are now ready to numerically compute the finite basis approximation to the 
Laplace operator (/a|A|/b) and the coefficient matrix {fa\fb) for any fixed values 
of and n^. The coefficients do not form the unit matrix, indicating that the 
approximating basis eq. fl30p of 3^k^ is not orthonormal. Even though one could 
orthonormalize the basis, this would be numerically unsound and it is easier to directly 
solve the generalized eigenvalue problem eq. f|TT]) . We implemented this algorithm in 
C++. In practice, the most time-consuming part is the evaluation of the numerical 
integrals for the matrix elements of the Laplace operator. We perform this step 
in parallel on a 10-node dual Opteron cluster, using MPI [IS] for communication. 
Finally, we use LAPACK [17] to compute the eigenvalues and eigenvectors. Note 
that the matrix eigenvectors are the coefficients {fa\4>) and, hence, the corresponding 
eigenfunction is 

dim3"fe^-l 
a=0 

We present our results in two ways. First fix A;^, thus restricting the total number of 
non-degenerate eigenvalues A„ to dimJ'fc^. These eigenvalues are then plotted against 
the number of points that we use to evaluate an integral. For smaller values of 
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Figure 1: Spectrum of the scalar Laplacian on with the rescaled Fubini- 
Study metric. Here we fix the space of functions by choosing degree 
= 3, and evaluate the Laplace operator at a varying number of 
points n^. 



n^, the eigenvalues are fairly spread out. However, as is increased the eigenvalues 
break into distinct groups, each of which rapidly coalesces toward a unique value. 
One can then compare the limiting value and multiplicity of each group against the 
exact analytic result. We find perfect agreement. To be concrete, let us present 



the numerical results for the case = 3. We plot these results in Figure 1 As 



is increased from 10,000 to 1,000,000, the dime's = 400 eigenvalues A„ cluster 
into 4 distinct groups with multiplicity 1, 15, 84 and 300. These clusters approach 
the theoretical values of the first four eigenvalues respectively, as expected. That 
is, the numerically calculated eigenvalues condense to the analytic results for the 
eigenvalues and multiplicities listed in lTableH on page [91 At any n^, the eigenfunction 
(pn associated with each A„ is evaluated as a sum over the basis functions {/a|a = 
0, . . . , 399}. We do not find it enlightening to present the numerical coefficients. 



13 



300 



250 



200 



150 



100 



50 



/is = 1911 




/i4 = 825 






/i3 = 300 










Ai2 = 84 








/ii = 15 







276.622 



193.635 



124.48 

69.1554 

27.6622 




1 2 3 4 5 

Figure 2: Spectrum of the scalar Laplacian on with the rescaled Fubini- 
Study metric. Here we evaluate the spectrum of the Laplace opera- 
tor as a function of k^, while keeping the number of points fixed at 
= 100,000. Note that k^ determines the dimension of the matrix 
approximation to the Laplace operator. 



The second way to present our numerical results is to fix and study the de- 
pendence of the eigenvalues on k^p. As was discussed in [Subsection 3.1t since the 
eigenfunctions of the Laplace operator are linear combinations of the elements of our 
basis, the accuracy of A„ should not depend on k^. However, increasing k^ does add 
higher-frequency functions to the approximating space of functions. More explicitly, 
going from k^p to /c^+i will add an extra fik^_^_^ eigenvalues to the numerical spec- 
trum, corresponding to the dimension of the (/c^+i, 0, fc^+i) irreducible representation 
of SU{4). This is exactly the behavior that we observe in [Figured 
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3.3 Asymptotic Behaviour 



It is of interest to compare the asymptotic behaviour of the numerical solution to 
the theoretical prediction of Weyl's formula, which determines the asymptotic growth 
of the spectrum of the scalar Laplace operator. Specifically, it asserts that on a 

2 

Riemannian manifold X of real dimension (i, the eigenvalues grow as A„ ~ ra^ for 
large n. Here it is important to keep track of multiplicities by including the degenerate 
eigenvalue multiple times in the sequence {\n\i as we do in our numerical calculations. 
The precise statement of Weyl's formula is then that 

xT (47r)fr(f + l) 

tL^= voi(x) • 

Applying this to P^, which has d = Q and the volume scaled to Volx(P'^) = 1, we find 
that 

\3 

lim = 3847r3. (37) 



In Figure 3 we choose /c,^ = 3 and plot ^ as a function of n for the numerical 



values of A^, as well as for the exact values listed in ITable II The numerical results 
are presented for six different values of n^. For each value of n^, as well as for the 
exact result, the — break into three groups, corresponding to the first three massive 
levels with multiplicities 15, 84, and 300, respectively. Note that, as gets larger, 
the numerical results converge to the exact result. That is, each segment approaches 
a curve of the form ^^^st. Furthermore, as the number of eigenvalues increase, the 
end-points of the curves asymptote toward the Weyl limit 3847r^. 



4 Quintic Calabi-Yau Threefolds 

Quintics are Calabi-Yau threefolds Q C F^. Denote the usual homogeneous coordi- 
nates on by z = [zq : Zi : Z2 '■ Z3 : z^]. A hypersurface in P^ is Calabi-Yau if and 
only if it is the zero locus of a degree-5 homogeneous polynomial 

Qiz) = C(no,ni,n2,n3,n4)^0°^1^^2^^3^^4''- (38) 

'^0+™i+"-2+n.3+"-4=5 

By the usual abuse of notation, we denote both the defining polynomial Q{z) and 
the corresponding hypersurface {Q{z) = 0} C P'^ by Q. There are (^'''4"^) = 126 
degree-5 monomials, leading to 126 coefficients C(^no,ni,n2,n3,n4,) ^ C. These are not all 
independent complex structure parameters, since the linear GL(5, C)-action on the 
five homogeneous coordinates is simply a choice of coordinates. Hence, the number 
of complex structure moduli of a generic quintic Q is 126 — 25 = 101. 
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Figure 3: Check of Weyl's formula for the spectrum of the scalar Laplacian 
on with the resettled Fuhini- Study metric. We fix the space of 
functions by taking = 3 ttnd evaluate ^ as a function of n at 



a varying number of points n^j,. 
eigenvalues is the same as for k^ 



Note that the data used for the 
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A natural choice of metric on P is the Fubini-Study metric gij = did^K-ps, where 



I 4 

i^Fs = - In \^ ZiZj. (39) 



'■FS 

TT 

i=0 



This induces a metric on the hypersurface Q, whose Kahler potential is simply the 
restriction. Unfortunately, the restriction of the Fubini-Study metric to the quintic 
is far from Ricci-flat. Recently, however, Donaldson [31] presented an algorithm for 
numerically approximating Calabi-Yau metrics to any desired accuracy. To do this 
in the quintic context, one takes a suitable generalization, that is, one containing 
many more free parameters, of the Fubini-Study metric. The parameters are then 
numerically adjusted so as to approach the Calabi-Yau metric. 

Explicitly, Donaldson's algorithm is the following. Pick a basis for the quotient 

C[z,,...,z,],/{Q{z)) (40) 

of the degree-/c polynomials on P"^ modulo the hypersurface equation. Let us denote 
this basis by Sq,, a = 0, . . . , N{k) — 1 where 

('+^^) < A; < 5 

Nik) = I ^ (41) 

For any given quintic polynomial Q{z) and degree /c, computing an explicit polynomial 
basis {sq} is straightforward. Now, make the following ansatz 

N{k)-1 
a,/3=0 

for the Kahler potential. The hermitian N{k) x A^(/c)-matrix parametrizes the 
metric on Q and is chosen to be the unique fixed point of the Donaldson T-operator 

N{k) f SaSj^ 



Tih)-3 = / ^ ^ dVolcY, (43) 

where 

dVolcY = A (44) 

and Q is the holomorphic volume form. The metric determined by the fixed point 
of the T-operator is called "balanced". Hence, we obtain for each integer k > 1 the 
balanced metric 

_ N{k)-1 
o,/3=0 



17 



Note that they are formally defined on but restrict directly to Q, by construction. 
One can show [IE] that this sequence 



(fc) 
3ij 



(46) 



of balanced metrics converges to the Calabi-Yau metric on Q. 

It is important to have a measure of how closely the balanced metric g\^'^ at a given 
value of k approximates the exact Calabi-Yau metric . One way to do this is the 

following. Let g^^^ be a balanced metric, Uk the associated (1, l)-form and denote by 



Vol^(Q,fc) = /J^, VolcY(g) 



1] Afi 



(47) 



the volume of Q evaluated with respect to Uk and the holomorphic volume form ^2 
respectively. Now note that the integral 



(rk{Q) 



1 



voicY (g) Jq 



^/VoW(Q,A;) 



QhQ VolcY [Q) 



dVoI 



CY 



(4^ 



must vanish as ujk approaches the Calabi-Yau Kahler form. That is 



0. 



(49) 



Following [33], we will use ak as the error measure for how far balanced metric g^~^ 
is from being Calabi-Yau. Finally, to implement our volume normalization we will 
always scale the balanced metric so that 



Vol,, {Q,k)=l 



(50) 



at each value of k. 



4.1 Non-Symmetric Quintic 

In this subsection, we will pick randonj^ coefficients C(^no,ni,n2,ni,n4) for the 126 different 
quintic monomials in the 5 homogeneous coordinates. An explicit example, which we 
use for the analysis in this section, is given by 

Q[z) = (-0.319235 + 0.709687i)^^ + (-0.327948 + 0.811936i)42i 

+ (0.242297 + 0.219818z)4^2 + ■ ■ ■ + (-0.265416 + 0.122292z)zt (51) 

^To be precise, we pick uniformly distributed random numbers on the unit disk {zgC:|z|<1}. 
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We refer to this as the "random quintic". Of course, any other random choice of 
coefficients would lead to similar conclusions. The polynomial eq. (15T1) completely 
fixes the complex structure. Furthermore, the single Kahler modulus determines the 
overall volume, which we set to unity. 

Using Donaldson's algorithm [311 [33], HQ] which we outlined above, one can compute 
an approximation to the Calabi-Yau metric on the quintic defined by eq. ( l5Tl) . The 
accuracy of this approximation is determined by 

• The degree k G Z>o of the homogeneous polynomials used in the ansatz eq. f|^2|) 
for the Kahler potential. To distinguish this degree from the one in the approx- 
imation to the Laplace operator, we denote them from now on by kh and /c<^, 
respectively. In this section, we will use 

kh = 8. (52) 

Note that the choice of degree kh determines the number of parameters 

/i"^ G Mat {N{kh) X N{kh), C) (53) 

in the ansatz for the Kahler potential, eq. (142|) . This is why kh is essentially 
limited by the available memory. We choose kh = 8 because it gives a good 
approximation to the Calabi-Yau metric, see below, without using a significant 
amount of computer memory (^ 7 MiB). 

• The number of points used to numerically integrate within Donaldson's T- 
operator [33]. To distinguish this number from the number of points used to 
evaluate the Laplacian, we denote them by and Uh respectively. As argued 
in [10], to obtain a good approximation to the Ricci-fiat metric one should choose 
Uh 3> N{khYi where N{kh) is the number of degree- /c/i homogeneous monomials 
in the 5 homogeneous coordinates modulo the Q{z) = constraint, see eq. (14T!) . 
In our computation, we will always take 

n,, = 10 ■ N{khf + 50,000. (54) 

This rather arbitrary number is chosen for the following reasons. First, the 
leading term assures that Uh ^ N{kh)^ by an order of magnitude and, second, 
the addition of 50,000 points guarantees that the integrals are well-approximated 
even for small values of kh- It follows from eq. (HTj) that for kh = 8 we will use 

Uh = 2,166,000 (55) 

points in evaluating the T-operator. 
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Using the Donaldson algorithm with kh and Uh given by eqns. (1521) and ( !55i) respec- 
tively, one can now compute a good approximation to the Calabi-Yau metric in a 
reasonable amount of tim^. The expression for the metric itself is given as a sum 
over monomials on Q of degree kh = 8 with numerically generated complex coeffi- 
cients. It is not enlightening to present it here. However, it is useful to compute the 
error measure defined in eq. ( I48l) for this metric. We find that 

cTg ^ 5 X 10"^ (56) 

meaning that, on average, the approximate volume form -gf and the exact Calabi-Yau 
volume form f2 A agree to about 5%. Finally, having found an approximation to 
the Ricci-fiat metric, one can insert it into eq. ([1]) to determine the form of the scalar 
Laplacian. 

We can now compute the spectrum of the scalar Laplace operator as discussed in 
the previous section. First, one must specify a finite-dimensional approximation to 
the space of complex- valued functions on Q. For any finite value of ks, we choose 



J'fc = span 



S(yS 



aJ = 0,...,Nik^)-l}, (57) 



where {sa\a = 0, . . . , N{k^) — 1} are a basis for the homogeneous polynomials modulo 
the hypersurface constraint 

span{s„} = C [zo, z^]^^ j {Q{z)). (58) 

Such a basis was already determined during the Donaldson algorithm for the metric, 
the only difference being that now the degree is k^ instead of k^. The counting function 
N{k^) is given by eq. (l4T|l . Clearly, 

c\im-Jk, = N{k^f. (59) 

Computing the matrix elements of the Laplace operator requires another numerical 
integration which is completely independent of the one in the T-operator. We denote 
the number of points in the matrix element integration by n^, as we did in the previous 
section. We first present the resulting eigenvalue spectrum for fixed fc^ = 3 plotted 



against an increasing number of points n^. Our results are shown in Figure 4 From 
eq. ( l4TI) we see that A^(3) = 35 and, hence, there are 35^ = 1,225 non-degenerate 
eigenvalues Aq, . . . , Ai,224- Note that for smaller values of the eigenvalues are 
fairly spread out, and that they remain so as Ucj, is increased. This refiects the fact 
that for any Calabi-Yau manifold there is no continuous isometry, as there was for 
the P^. Furthermore, for the random quintic eq. (ISTj) there is no finite isometry 



'^That is, within a few hours of "wall" time. 
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Figure 4- Eigenvalues of the scalar Laplace operator on the same "random 
quintic" defined in eq. (j5ip . The metric is computed at degree kh = 8, 
using rih = 2,166,000 points. The Laplace operator is evaluated at 
degree = 3 on a varying number of points. 
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group either. Therefore, one expects each eigenvalue to be non-degenerate, and our 
numerical results are clearly consistent with this. At any n^, the eigenf unctions 0„ 
are a linear combination of the 1,225 basis functions. We do not find it enlightening 
to list the numerical coefficients explicitly. 

Note that the accuracy of the numerical integration for the matrix element^ is 
not as crucial as in the T-operator, since we are primarily interested in the low lying 
eigenvalues corresponding to slowly-varying eigenfunctions. This is nicely illustrated 
by Figure 4 where the eigenvalues rather quickly approach a constant value as we 
increase n^, even though <^ rih- For this reason, = 200,000 gives a sufficiently 
good approximation and we will use this value for the reminder of this subsection. 

A second way to present our numerical results is to fix ra^ and study the dependence 
of the eigenvalues on /c^. This is presented in Figure 5 We first note that the number 
of eigenvalues indeed grows as Nlk^)"^, as it must. Second, as one expects, the smaller 
eigenvalues do not change much as one increases fc^. The higher eigenvalues, however, 
depend strongly on the truncation of the space of functions, since their eigenfunctions 
vary quickly. 

Finally, we plot A^/n as a function of n in Figure 6 We see that this ratio does 
approach the theoretical value of 3847r'^ as and n increase. This confirms that 
the volume normalization in eq. (150|) is being correctly implemented and that our 
numerical results are consistent with Weyl's formula eq. (l36l) . 



4.2 Fermat Quint ic 

We repeat the analysis of the previous section for the Fermat quintic defined by 

Qf{z) = + zl + zI + zI + z\. (60) 

As before, the single Kahler modulus is chosen so that the volume of the Fermat 
quintic is unity. Now, however, we are at a different point in the complex structure 
moduli space, eq. fl^Ul) instead of the random quintic eq. fl^T]) . Hence, we will perform 
the numerical integrations now using points lying on a different hypersurface inside 
P^. Except for using different points, we compute the Calabi-Yau metric ov^Qp using 
Donaldson's algorithm exactly as in the previous subsection. In particular 

• The degree G Z>o of the homogeneous polynomials used in the ansatz eq. (H2|) 
for the Kahler potential is chosen to be 

kh = 8. (61) 

This is the same degree as we used for the random quintic. 

^Recall that — + oo is the continuum Umit for the numerical integration in the T-operator, and 
— > oo is the continuum hmit for the numerical integration determining the matrix elements of the 
Laplace operator. 
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Figure 5: Eigenvalues of the scalar Laplace operator on a random quintic plot- 
ted against k^. The metric is computed at degree = 8, using 
rih = 2,166,000 points. The Laplace operator is then evaluated at 
= 200,000 points. 
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Figure 6: Check of Weyl 's formula for the spectrum of the scalar Laplace oper- 
ator on a random quintic. The metric is computed at degree kh = 8, 
using = 2,166,000 points. The Laplace operator is evaluated at 
riff, = 200,000 points and degrees k^i) = 1,2,3. Note that the data 
for the eigenvalues is the same as in Figure 5[ According to Weyl's 



formula, the exact eigenvalues have to satisfy lim A„/n = 3847r 
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• We take the number of points used to numerically integrate Donaldson's T- 
operator to be 

n;, = 10 ■ N{8Y + 50,000 = 2,166,000 (62) 

This satisfies the condition that ^ NikhY, ensuring that the numerical 
integration is sufficiently accurate. 

Using kh and rih given by eqns. (16T|) and (162|) respectively, one can compute an ap- 
proximation to the Calabi-Yau metric using Donaldson's algorithm. The numerical 
expression for the metric is tedious and will not be presented here. The error measure 
eq. (HHl) for this kh = 8 balanced metric is 

as ^ 5 X 10"^ (63) 

Hence, the approximate volume form gf and the exact Calabi-Yau volume form Afi 
agree to about 5%. The metric determines the scalar Laplacian, eq. ([T]). 

To determine the matrix elements of the Laplace operator, one has to select an 
approximating basis for the linear space of complex functions on Qp, eq. fl60|) . For any 
finite /c0, we again choose the function space 5"/;^ as in eqns. fl57j) and fl58l) . This basis 
was already determined during the Donaldson algorithm for the metric. Computing 
the matrix elements of the Laplace operator requires another numerical integration 
which is completely independent of the one in the T-operator. As we did previously, 
we denote the number of points in the matrix element integration by n^. 

We first present the resulting eigenvalue spectrum for fixed k^ = 3 plotted against 



an increasing number of points n^. Our results are shown in Figure 7 Note from 
eq. 0591) that the total number of eigenvalues is given by dime's = A^(3)^ = 1,225. 
One immediately notices a striking difference compared to the analogous graph for 



the random quintic, [Figure 4} Here, the eigenvalues converge towards degenerate 
levels. For smaller values of n<^, the eigenvalues are fairly spread out. However, as 
n<^ is increased the eigenvalues begin to condense into degenerate levels. Clearly, this 
must be due to symmetries of the Fermat quintic. As mentioned above, no Calabi- 
Yau manifold has a continuous isometry. However, unlike the random quintic, the 
Fermat quintic eq. (!60|l does possess a finite isometry group, which we will specify 
below in detail. Therefore, the exact eigenvalues of A on should be degenerate 
with multiplicities given by the irreducible representations of this finite group. As we 
will see in [Subsection 4.31 the numerically computed degeneracies of the eigenvalues 
exactly match the irreducible representations of a this finite isometry group. Again, 
we do not find it enlightening to present the numerical results for the eigenfunctions. 
Moreover, as discussed previously, the accuracy of the matrix element integration for 
low-lying eigenvalues need not be as great as for the T-operator. As is evident from 



Figure 7\ a value of = 500,000 is already highly accurate and we will use this value 



in the remainder of this subsection. 
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Figure 7: Eigenvalues of the scalar Laplace operator on the Fermat quintic. 

The metric is computed at degree = 8, using Uf^ = 2,166,000 
points. The Laplace operator is evaluated at degree = 3 using a 
varying number of points. 



26 



200 



^ 150 



100 



50 







= 30 


/i4 


= 60 




= 4 


/i2 


= 20 



fii = 20 



Figure 8: Eigenvalues of the scalar Laplace operator on the Fermat quintic. 

The metric is computed at degree = 8, using = 2,166,000 
points. The Laplace operator is evaluated at = 500,000 points 
with varying degrees k^. 
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A second way to present our numerical results is to fix as in the previous 
paragraph and study the dependence of the eigenvalues on /c^. This is presented in 
Figure 8 We first note that the number of eigenvalues grows as N{k^Y, as it must. 



Second, as one expects, the smaller eigenvalues do not change much as one increases 
/c<^, whereas the higher eigenvalues depend strongly on the truncation of the space of 
functions. This is also to be expected, since their eigenf unctions vary quickly. 
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Figure 9: Check of WeyVs formula for the spectrum of the scalar Laplace oper- 
ator on the Fermat quintic. The metric is computed at degree = 8, 
using Uh = 2,166,000 points. The Laplace operator is evaluated at 
riff) = 500,000 points and degrees = 1,2,3. Note that the data 



for the eigenvalues is the same as in Figure 8. According to WeyVs 
formula, the exact eigenvalues have to satisfy lim A^/n = 3847r^. 



Third, let us plot A^/n as a function of n in Figure 9 
the theoretical value of 3847r'^ 



This ratio approaches 



as k^ and n increase. This confirms that the volume 
normalization in eq. fl5U]l is being correctly implemented and that our numerical results 
are consistent with Weyl's formula eq. 



28 



4.3 Symmetry Considerations 

Recall from [Figure 7\ that the eigenvalues of the scalar Laplace operator condense 
to a smaller number of degenerate levels as — cxo, that is, in the limit where 
the numerical integration becomes exact. The same phenomenon is clearly visible 



at different values of k^i,, see [Figure 8| Of course the eigenvalues are never exactly 
degenerate due to numerical errors, but counting the nearby eigenvalues allows one to 
determine the multiplicities. Averaging over the eigenvalues in each cluster yields an 



approximation to the associated degenerate eigenvalue. Using the data from Figure 8^ 
we list the low-lying degenerate eigenvalues and their multiplicities in ITable 2[ As 



m 





1 


2 


3 


4 


5 




1.18 X 10"^^ 


41.1 ±0.4 


78.1 ±0.5 


82.1 ±0.3 


94.5 ± 1 


102 ± 1 




1 


20 


20 


4 


60 


30 



Table 2: The degenerate eigenvalues Am and their multiplicities fim on the Fer- 
mat quintic, as computed from the numerical values calculated with 
kh = 8, n/i = 2,166,000, = 3, = 500,000. The errors are the 
standard deviation within the cluster of fin numerical eigenvalues. 

discussed previously, multiplicities in the spectrum of the Laplace-Beltrami operator 
results must follow from some symmetry. In [Section 31 we saw that the SU (4) symme- 
try of leads to degenerate eigenspaces of the scalar Laplacian. However, a proper 
Calabi-Yau threefold never has continuous isometries, unlike projective space. Never- 
theless, a suitable non-AbeliaiJ^ finite group action is possible and, in fact, explains 
the observed multiplicities, as we now show. 

First, note that for each distinct eigenvalue the corresponding space of eigenfunc- 
tions must form a representatior0 of the symmetry group. Clearly, the degeneracies 
of the eigenvalues observed in [Figure 7" and Figure 8 must arise from an isometry of 



Qp. In fact, the Fermat quintic does have a large non-Abelian finite symmetry group. 
To see this, note that the zero set of eq. (!60|) is invariant under 

• Multiplying a homogeneous coordinate by a fifth root of unity. However, not 
all (Zs)^ phases act effectively because the projective coordinates are identified 



^Interestingly, the correct multiplicity /ii = 20 was derived by a completely different argument 
in [49]. 

"'^'^An Abelian symmetry group would only have one-dimensional representations and, hence, need 
not lead to degenerate eigenvalues. Note that any finite group has a finite number of irreducible 
representations and, therefore, one expects only a finite number of possible multiplicities for the 
eigenvalues of the Laplace operator. This is in contrast to the aforementioned P"^ case, where the 
multiplicities grow without bound. 

^^An actual linear representation, not just a representation up to phases (projective representa- 
tion). 
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d 


1 2 


4 


5 6 8 


10 


12 20 30 40 60 


80 120 


# of irreps 
in dim d 


4 4 


4 


4 2 4 


4 


2 8 8 12 18 


4 2 



Table 3: Number of irreducible representations o/Aut((5i?) = Z2 x Aut((5i;') 
in eac/i complex dimension. 

under the rescaling 

[zq : 21 : : ^3 : ^4] = [A^q : A^i : \z2 : A^s : Xz^ . (64) 
Only {l^f/Z^ ~ (Zg)^ acts effectively. 

• Any permutation of the 5 homogeneous coordinates. The symmetric group 6*5 
acts effectively. 

• Complex conjugation Z2. 

The first two groups act by analytic maps, and together generate the semidirect 
product 

Aut [Qf) = ^5 k (Zs)^ (65) 
of order 75, 000. Our notation and the relevant group theory is discussed in Appendix B 



The full discrete symmetry group, including the complex conjugation Z2, is 

^^{Qf) = Z2 k Aut [Qf] = (55 X Z2) k {Z^f (66) 

and of order 150,000. Note that even though the Z2 acts as complex conjugation 
on the base space, the whole Ant{QF) acts linearly on the the basis of complex 
functions on Qf and, hence, on the eigenf unctions. There are 80 distinct irreducible 
representations occurring in 14 different dimensions, ranging from 1 to 120. We list 
them in ITable 31 

We conclude by noting that the multiplicities listed in ITable 21 also occur in 
ITable 3[ That is, the eigenspaces of the degenerate eigenvalues of the scalar Laplacian 
on Qf, computed using our numerical algorithm, indeed fall into irreducible repre- 
sentations of the finite symmetry group (5*5 x Z2) x (Z5) , as they must. This gives 
us further confidence that our numerical computation of the Laplacian spectrum is 
reliable. 



4.4 Donaldson's Method 

Donaldson [31] conjectured a method to compute the eigenvalues of the scalar Laplace 
operator that is completely independent of our approach. His calculation of the spec- 
trum of the scalar Laplacian is very much tied into his algorithm for computing 
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balanced (Calabi-Yau) metrics. In our algorithm, on the other hand side, any metric 
could be used and no particular simplifications arise just because the metric happens 
to be balanced or Calabi-Yau. Because they are so different, it is quite interesting 
to compare both methods. We will now review his proposal, and then compare it 
with our previous computation of the eigenvalues on the Fermat quintic as well as the 
random quintic. 

In this alternative approach to calculating the spectrum of the Laplace-Beltrami 
operator, one first has to run through Donaldson's algorithm for the metric. In par- 
ticular, one had to choose a degree k, fix a basis {s^la = 0, . . . , N{k) — 1}, and obtain 
the balanced metric /i"^ as the fixed point of Donaldson's T-operator. Let us write 

(^"'M = V^#- (67) 



for the integrand of the T-operator eq. P3l) . Donaldson's alternative calculation of 
the eigenvalues then hinges on the evaluation of the integral 



X 



where we again normalize Vol(X) = 1. One can think of Q as a linear operator on 
the space of function^ 

5-° = span|(s„,s^) 0<aJ<N{k)-iy (69) 



acting via 

Q : 3^° ^ (s^.sp) ^ Y.Q--Msh^''h'\s,,s^). (70) 
In [21], Donaldson conjectures that 



lim Q = e ^-Ww (71) 

fc— >oo 

as operators on 

lim J° = C°^(X,C). (72) 

fc— >oo 

For explicitness, let us look in more detail at the individual steps as they apply to 
any quintic X = Q C P^: 

1. First, pick a degree k and a basis [^sq, . . . , SN{k)-i] of degree-/c homogeneous 
polynomials modulo the hypersurface equation Q = Q. 



^^Note the similarity with the approximate space of functions 'Jk^ used previously, eq. ([57)) . When 
computing the matrix elements of the Laplace operator directly, the precise form of the denominator 
is not overly important as long as it has the correct homogeneous degree, and we always chose 
(S NiP)*^* foi' simplicity. 
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2. Compute the Calabi-Yau metric via Donaldson's algorithm. It is determined by 
the N{k) X N{k) hermitian matrix h°'^. 

3. Compute the N[kY scalar integrals in eq. (|68l) . The numerical integration can 
be performed just as in Donaldson's T-operator, see ISection 4l 

4. Compute the N{kY x N{kY matrix 

N{kh)-l 

QNik).+f^'^'^' = E Qo^-p^^rh^'h-^' (73) 

CT,r=0 

and find its eigenvalues A„. Note that Q{ is not hermitiaiJ^ and one should use 
the Schur factorizatioE0 to compute eigenvalues. 

5. Discard all A„ < 0, these correspond to high eigenvalues of the Laplacian that 
are not approximated well at the chosen degree k. The eigenvalues of the scalar 
Laplace operator are 

A„ = -ATi^/Nik) lnA„. (74) 

We note that, in this approach to the spectrum of the Laplace-Beltrami operator, 
there is only one degree k that controls the accuracy of the eigenvalues of the scalar 
Laplacian and at the same time the accuracy of the Calabi-Yau metric. In fact, com- 
puting the integral eq. at degree k is about as expensive as computing Donaldson's 
T-operator at degree 2k. In other words, a general limitation of this approach is that 
one has to work with a relatively low precision metric. 



In Figure 10 we compare the two approaches for computing the spectrum of the 
Laplace-Beltrami operator on the Fermat quintic. We compute the eigenvalues using 
Donaldson's method at degrees k = 1,2,3 and evaluate the necessary integral eq. ( l68l) 
using n = lON^k)'^ + 100,000 points. For comparison, we also plot the eigenvalues 
obtained by directly computing the matrix elements of the Laplacian which we always 
compute at degree k^ = 3 using = 500,000 points. To estimate the effect of 
the metric on the eig envalues, we run our algorithm first with the metric obtain at 
degree kh = S anqlj Uh = 62250 (bad approximation to the Calabi-Yau metric, red 
diagonal crosses) as well as with kh = S and Uh = 2,166,000 (good approximation 
to the Calabi-Yau metric, blue upright crosses). We find that the eigenvalues do not 
strongly depend on the details of the metric. Generally, Donaldson's method and the 
direct computation yield very similar results. There is a slight disagreement for the 
second and third massive level, where the matrix element calculation points toward 
/i2 = 20, /i3 = 4 while Donaldson's method suggests the opposite order /is = 4, /i4 = 
20. We suspect this is to be a numerical error due to the finite degrees and it would 
be interesting to go to higher degree in k, ks, and k^,. 



^^Ql is, however, conjugate to a hermitian matrix and hence has real eigenvalues. 
"'^^Instead of the dqds algorithm we use for computing eigenvalues of hermitian matrices. 
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The number of points nu is always obtained from the heuristic eq. ([5^ . 
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Figure 10: Donaldson's method of computing the spectrum (polygon symbols) 
of the scalar Laplace operator on the Fermat quintic compared to 
our direct computation (crosses). Note that the blue symbols are 
the highest- accuracy values, respectively. See Subsection 4-4\ for 
further discussion. 
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Figure 11: Donaldson's method of computing the spectrum (polygon symbols) 
of the scalar Laplace operator on the random quintic compared to 
our direct computation (crosses). Note that the blue symbols are 
the highest- accuracy values, respectively. In Donaldson's method 
the numerical integration was performed with n = 10A^(A;) + 100,000 
points. In the direct computation, the metric was approximated at 
degree = 8 using = 2,166,000 points and the Laplace operator 
was evaluated at = 500,000 points. 
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Finally, in [Figure 11] we repeat this comparison for the quintic eq. (1511) with ran- 
dom coefficients. In this case, there are no discrete symmetries and one expects all 
massive levels to be non-degenerate. We again find good agreement between the two 
approaches towards solving the Laplace equation. 



5 Z5 X Z5 Quotients of Quintics 

Thus far, we have restricted our examples to quintic Calabi-Yau threefolds Q C P'*. 
These manifolds are simply connected by construction. However, for a wide range 
of applications in heterotic string theory we are particularly interested in non-simply 
connected manifolds where one can reduce the number of quark/lepton generations 
as well as turn on discrete Wilson lines. Therefore, in this section we will consider 
the free Z5 x Z5 quotient of quintic threefolds, see [40] for more details. 

5.1 Z5 X Z5 Symmetric Quintics and their Metrics 

Explicitly, the group action on the homogeneous coordinates [zq : ■ ■ ■ : z^^ E is 

gi : [zo ■■ Zi : Z2 : Z3 : z^l — > [zq : e ^ Zi : e ^ Z2 : e z^ : e ^ zA, 

1 ^^^^ 
92 ■ \Zq . Zi . Z2 ■ z^ . 2:4] — > yzi .Z2.Z3.Z4. zqJ • 

As we discussed in [Section 4l a generic quintic is a zero locus of a degree-5 polynomial 
containing 126 complex coefficients. However, only a small subset of these quintics is 
invariant under the Z5 x Z5 action above. As we will show below, the dimension of the 
space of invariant homogeneous degree-5 polynomials is 6. Taking into account that 
one can always multiply the defining equation by a constant, there are 5 independent 
parameters (pi, ... (p^ & C Thus, the Z5 x Z5 symmetric quintics form a five parameter 
family which can be written as 

Q{z) = {zl + z\ + zl + zl + ^1) 

+ (\)x{zqZiZ2Z3Z^ 

( 3 3 3 3 3 \ 

ZqZiZ4 + ZqZ^Z2 + ZqZ^Z^ + Z1Z2Z3 + Z2Z3Z4) 

,1/2 2,2 2,2 2,2 2,2 2\ \'^) 
-f- (p3yZQZiZ2 -r Z1Z2Z3 + Z2Z3Z^ -\- Z3Z4ZQ + Z^ZqZiJ 

,1/22 ,2 2 ,2 2 ,2 2 , 2 2 \ 
-T (PiyZ^Z-^Z^, -f- Z^Z2Z4 -f- Z2Z3ZQ -\- Z^Z^Zi -\- Z^ZqZ2J 

+ 05(2^0^2^3 + zfz^Z^ + Z2Z4ZQ + Z^ZqZi + zlziZ2), 

where 0i, . . . , 05 G C are local coordinates on the complex structure moduli space. 
From now on, Q C P^ will always refer to a quintic of this form. 
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For generic coefficient^^ 0j, the hypersurface Q is a smooth Calabi-Yau threefold. 
Moreover, although the group action eq. (1751) necessarily has fixed points in P^, these 
fixed points do not intersect a generic hypersurface Q. Thus the quotient 

g = g/(Z5xZ5) (77) 

is again a smooth Calabi-Yau threefold. As a general principle, we will compute 
quantities on the quotient Q by computing the corresponding invariant quantities on 
the covering space Q. For example, the complex structure moduli space of Q is the 
moduli space of Z5 x Zs-invariant complex structures on Q. Hence, its dimension is 

/i2'i(g) = ^mvH^'\Q) = dimi72.i(g)Z5xZ5 ^ 5^ ^^g) 

corresponding to the 5 independent parameters 0i, . . . , 05 in a Z5 x Zs-invariant 
quintic Q{z). 

In the same spirit, we will compute the Calabi-Yau metric on Q by performing the 
analogous computation on the covering space Q. To begin, one must choose a degree 
kh and determine a basis Sa for the corresponding Z5 x Zs-invariant homogeneous 
degree- fcft polynomials 

span{. J = C[^o, • • • , ^4]fr^y (Q(^)> (79) 
on Q. Note, however, that for any homogeneous degree-Zc^ polynomial Pk^z) 

91929,' 9^\pk,X^)) = e'^^^ PkM (80) 

and, hence, the two Z5 generators in eq. fl75|) do not always commute. It follows that 
for a space of homogeneous polynomials to carry a linear representation of Z5 x Z5, 
let alone have an invariant subspace, their degree kh must be divisible by 5; that is, 

kh e 5Z. (81) 

This can be understood in various ways, and we refer to [ID] for more details. Hence- 
forth, we will assume that eq. fIST]) is satisfied. 

The ffist step in determining the basis of sections {sq} on Q is to find a basis 
for the invariant polynomials C[zo, . . . ,^4]^^^^^ on P^. Such a basis is given by the 
Hironaka decomposition 

100 

C[Zo, ^1, Z2, Z3, Z^fkl"^' = r/, C[^i, 62, ^3, ^4, ^5]fc,-deg(,,)- (§2) 

i=l 

^^For example, any sufRciently small neighbourhood of (0i, . . . , (p^) = (0, . . . , 0) £ C^. Note that 
setting all 4>i = yields the Fermat quintic Qf, see eq. ([60]) . 
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Here, the 9j = 6j{z) and rji = rii[z) are themselves homogeneous polynomials of 
various degree The 6'i, 9^ are called "primary invariants" and the rji, 
?7ioo are called "secondary invariants" . The primary and secondary invariants are not 
unique, but one minimal choice is [ID] 



and 



01 = z^ + zl + Z2+zl + z\ 

62 = ZoZiZ2Z^Z4 

03 = 20^1^4 + z\z2Zq + z'^Z'jZi + Z^Z4^Z2 + Z^ZqZ^ (§3) 

a _ ^10 I ^10 , ^10 I ^10 I ^10 
U4 — Zq -T Z-^ ~r ^2 + ^3 + 

/I 8 8 8 8 8 

C'5 = ZQZ2Z3 + Z1Z3Z4 + ^2^420 + Z^ZqZi + Z^ZiZ2 



Vi = 1, 

_2 2,2 2,2 2,2 2,2 2 
V2 — ZqZiZ2 + ZiZ2Z^ + ^2^3^4 + ^3 ^4^0 + ^A^OZi, 

_22 ,22 ,22 ,22 ,22 
?73 — ZqZ-^Z^ + 2;^ 2:2^4 + ^2^3^0 + ^3^4^! + Z^ZqZ2^ 

rji = zlz2Z3 + zfz^Zi + z\z^Zq + 2:32:0^1 + z\zxZ2, (84) 

''75 = ZqZ2 + 2]^ 23 +2:2^4 "I" ^3^0 "I" ^4^1? 



^7100 — + 2^;^ +2:2 +23 + 2:4 . 

For example, the 6-dimensional space of invariant degree-5 homogeneous polynomials 
on is 

100 

C[2;o,2;i,2;2,2;3,Z4]5'^^^ = C[^i, ^2, ^3, ^4, ^sls-dcgC,,) 

= r/i^iC © 77i^2C © r/i^3C © r]2C © r/3C © r/4C, 

thus proving eq. (!76l) . 

Using the Hironaka decomposition, we can now determining the basis in eq. (!79|) 
by modding out the equation Q{z) = which defines the covering space. This was 
discussed in |10] . The result is that one can simply eliminate the first primary invariant 
using 

Ol = -(pl92 - (P2O3 - 03^72 - 04?73 - 05^74, (86) 

yielding 

100 

span{sa} = ^Vi ^[6*2, 6*3, 6*4, 6'5]fch-deg(,?,) (87) 

i=l 



17 



The degrees of the 6j, rji are muhiples of 5, of course. 
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where a = 0, . . . ,N^^^^^{kh) - 1. The number N^'^^^^{kh) of Z5 x Zs-invariant ho- 
mogeneous degree-kh polynomials modulo Q = was tabulated in [ID]. In particular, 
the first three values are 

N^^^'^'^^O) = 1^ ^Z5xZ5(^5) ^ 5^ N^^x^^i^lO) = 35, (88) 

which we will use below. 

We now have everything in place to compute the metric on Q. First, one specifies 
the five complex structure parameters (pi which define the Z5 x Zs-symmetric covering 
space Q. Then, all one has to do is to replace the homogeneous polynomials in 
the procedure outlined in [Section 4l by Z5 x Zs-invariant homogeneous polynomials. 
Donaldson's algorithm then calculates the Calabi-Yau metric on the Z5 x Zs-symmetric 
quintic Q and, hence, the metric on the quotient Q = Q/ (Z5 x Z5). In fact, we use 
a refinement of this method which is even more efficient, that is, achieves higher 
numerical accuracy in less computing time. As it is not relevant to the spectrum 



of the Laplace operator, we relegate the details to Appendix C Henceforth, we will 



always use the following parameters in the computation of the metric. 

• The degree of the invariant homogeneous polynomials for the Kahler potential 
is taken to be 

kh = 10. (89) 

• The number of points used to evaluate the T-operator is 

n;, = 10 X of independent entries in h"^^ + 100,000 = 406,250. (90) 
Note that h"^^, the matrix of free parameters in Donaldson's ansatz for the met- 



ric, is block diagonal in Appendix C Therefore, the total number of independent 



entries is in fact 30,625 and not simply N^^^'^^ilOf = 1,225. 

As always, it is unenlightening to present the numerical result for the approximation 
to the Calabi-Yau metric. It is useful, however, to consider the error measure aiQ. As 
an important example, let us choose as our Calabi-Yau manifold the Z5 x Z5 quotient 
of the Fermat quintic Qp. The computation of the metric takes about half an hour 
of wall time, with the resulting error measure of aio = 2.8 x 10^^. 



5.2 The Laplacian on the Quotient 

Having computed the Calabi-Yau metric on the quotient Q = Q f (Z^ x Z5), we now 
turn to the calculation of the spectrum of the Laplace-Beltrami operator A. To begin, 
one must specify a finite-dimensional approximation to the space of complex valued 
functions on Q. Note, however, that the scalar functions on Q are precisely the 
invariant functions on the covering space Q. More formally, an invariant function on 
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Q is of the form q*f = foq, where / is a function on the quotient Q and q : Q ^ Q 
is the quotient map. Hence, we will specify a finite-dimensional approximation to the 
space of complex-valued Z5 x Zs-invariant functions on Q. For any finite value of k^, 
we choose 



SfyS 



«,^ = 0,...,iV^^^^^(A;^)-lL (91) 



where {sq,} is a basis for the invariant homogeneous polynomials modulo the hyper- 
surface constraint 

span{s J = C [zo, • • • , ^4]5'^-^ / {Q{z)). (92) 

We already had to determine such a basis while applying Donaldson's algorithm for 
the metric, the only difference now is that the degree is instead of k^. The counting 
function N^^^^^{k^) is the same, and some of its values were given in eq. (IHHjl . Clearly, 

dimlfg^^s = (^Ar^^^^5(A;^))'. (93) 

Having specified 3"^^^^^, we can now calculate any matrix element on Q simply by 

replacing the approximating space of functions on Q by the invariant functions on Q 
and integrating over Q. For example, the matrix elements of the Laplacian on Q are 

■„|A|/,> = / /,A/,dVol(g) = . ^ ^ I / (gVJA(g*/.)dVol(g). (94) 



Computing the matrix elements requires another numerical integration that is com- 
pletely independent of the one in the T-operator. As previously, we denote the number 
of points in the matrix element integration by n^. 

Having evaluated the matrix elements, one can now numerically solve the matrix 
eigenvalue equation eq. fITT]) for the eigenvalues and eigenfunctions of the Laplacian. 
Note that the factors of i^^^^^i cancel out of this equation, leaving identical eigenvalues 

and eigenfunctions on Q and Q, respectively. Since the functions in S^f^^^^ live on 
the covering space, we are actually solving 

AQ0n = A^^^^^0„, 0„ e C~(g, C)^^^^^ (95) 

on Q. Note that, as always, the volume measure of the integrals is chosen so that 
Vol(g) = 1. For the reasons stated above, the invariant eigenfunctions on Q can be 
identified with the eigenfunctions of the Laplacian on Q at the same eigenvalue, but 
with Vol(g) = i^^^^^i = ^. However, since we want to adhere to our convention 
of normalizing Vol(g) = 1, we have to rescale the volume and hence the eigenvalues 
Using eqns. ([3]) and (jll), the eigenvalues A„ on Q are 

\Z5xZ5 

(96) 



n 
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Using this method, one can compute the spectrum of the Laplace-Beltrami operator 
on the quotient of any Z5 x Z5 symmetric quintic. 



5.3 Quotient of the Fermat Quintic 

As an exphcit example, let us consider the quotient of the Fermat quintic. 
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Figure 12: Eigenvalues 
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of the scalar Laplace operator on the Fermat 
quintic Qp acting on Z5 x Z^-invariant eigenf unctions. Up to an 
overall factor due to our volume normalization, these are the same 
as the eigenvalues A.„ of the scalar Laplace operator on the quotient 
Qf = Qp/i'^b X ^5)- The metric is computed at degree = 10 
and niT^ = 406,250 points. The Laplace operator is evaluated using 
Ufj, = 100,000 points. 
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34.45 










2 
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101.2 
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Ai 


= 34.8 ±0.5 
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4 


4 


103.9 


35.53 










5 
6 


107.8 
109.2 


36.86 
37.36 


A2 


= 37.1 ±0.4 


/i2 = 


2 


7 


140.50 


48.05 










8 
9 


141.16 
141.47 


48.28 
48.38 


A3 


= 48.3 ±0.2 


/^3 = 


4 


10 


141.78 


48.49 










11 


149.57 


51.15 




A4 = 51.2 


yU4 = 


1 


12 
13 


"1 /^/^ n "1 

166.91 
169.48 


57.08 
57.96 


As 


= 57.5 ±0.6 


fJ'5 = 


2 


14 
15 


181.00 
184.15 


61.90 
62.98 


Ae 


= 62.4 ±0.8 


fJ'd = 


2 


16 


191.49 


65.48 










17 


193.55 


66.19 










18 


198.65 


67.94 











Table 4- Low-lying eigenvalues of the scalar Laplace operator on Qp, the Z5 x 
"L^-quotient of the Fermat quintic, computed with kh = k^f, = 10, 
Uh = 406,250, Utf) = 100,000. The first two columns are the numerical 
results. The third column specifies A, the average over the eigenvalues 
that are converging to a single degenerate level. The final column 
counts the multiplicities of each such level. 
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We numerically computed the spectrum of the scalar Laplace operator for each of 
the three values = 0,5,10 using eq. ( IHHl) . The resulting eigenvalues are shown 
Figure 12, Note that we present both the eigenvalues A^^^^^ on Q as well as the 



m 



normalized eigenvalues A„ on Q defined by eq. ( l96l) . 

We list the numerical values of the first few eigenvalues in ITable 41 and make the 
following two observations. First, the lowest eigenvalue Aq is no longer zero up to 
machine precision, as it was in ITable 2[ This is so because the constant function 
is not part of the approximate space of functions at /c^ = 10 and, therefore, the 
lowest eigenvalue Aq only approaches zero as increases. The actual numerical 
value Ao ~ 1.2 gives us an estimate of the error introduced by truncating the space 
of functions. Second, the low-lying eigenvalues clearly form degenerate levels. As 
usual, the numerical error caused by the truncation of the space of functions increases 
as we go to higher eigenvalues. However, the first 16 eigenvalues are sufficiently 
well separated that we can conjecture the underlying multiplicities /i. We list these 
degeneracies together with the best approximation to the true eigenvalue A in ITable ~^ 
Clearly, the degeneracies in the spectrum strongly hint at an underlying symmetry. 
We will discuss the associated isometry group in the following subsection. 

5.4 Group Theory and the Quotient Eigenmodes 

The free Z5 x Z5 action eq. (1751) is a subgroup of the symmetries of the Fermat quintic, 

Z5 X Z5 C AHt (Qi.) , (98) 

given in eq. fl66|) . Naively, one now would like to form the quotient to obtain the 
remaining symmetries on Qp = Qp/{Z^ x Z5). However, the Z5 x Z5 subgroup is 
not normal, that is, not closed under conjugation. The only possibility is to form the 
normal closurd^ 



(Z5 X Zs)^""*^^^^ = ^h^^gh ^ e Z5 X Z5, he Aut(gi.)}. (99) 

The quotient by the normal closure is well-defined, and we obtain 

ME{Qf) = AHt {Qf) I (Z5 X Zs)^"'^^^^ = Ao, (100) 

the dihedral group with 20 elements. However, just looking at the representation 
theory of Aut(Qir) is insufficient to understand the multiplicities of the eigenvalues 
of the Laplacian. Instead, one must use all of Aut(Qi;'), even those elements (called 
"pseudo-symmetries" in [H]) that do not correspond to symmetries of the quotient 
Qp. On a practical level, we also note that D20 has only 1- and 2-dimensional ir- 
reducible representations and could never explain the multiplicity [i\{Qp) = 4, for 
example, listed in ITable 4l 
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Also called the conjugate closure. 
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Table 5: Number nd of distinct irreducible representations ofAut{QF) in com- 
plex dimension d. We also list the dimension dim^^^^^ of the Z5 XZ5- 
invariant subspace for each representation, see eq. pU3|) . Note that it 
turns out to only depend on the dimension d of the representation. 



Qi 



Qi 



fiiiQp) = 20 

fi2{QF) = 20 

f^-siQp) = 4 
/X4(Qf) = 60 
fi^iQp) = 30 



fj'oiQp) 




f^2{Qp) 



4 
2 



Table 6: Projection of the multiplicity of eigenvalues on the Fermat quintic Q p 
to the Z5 X Z^-quotient Qp. 



As we discussed in ISubsection 4.3t the symmetry group of the Fermat quintic has 
80 distinct irreducible representations occurring in 14 different dimensions. Let us 
label them by pd,i, where d is the complex dimension and i = 1, . . . ,nd distinguishes 
the Ud different representations in dimension d. Under the Z5 x Z5 quotient 

Qf^Qp = Qp/i^6 X Z5) (101) 

all non-invariant eigenfunctions of the Laplacian are projected out and each invariant 
eigenfunction descends to an eigenf unction on Qp. Hence, the degeneracies of the 
eigenvalues are counted by the dimension 

dim (p^f (102) 



of the Z5 X Zs-invariant subspace. It turns out that, for the chosen Z5 x Z5 C Aut(Qi?), 
this dimension depends only on d, and not on the index i. We denote the common 
value by 

dim^^^^ = dim (pj-^') = . . . = dim (pj-/^) (103) 

and tabulate it in ITable 5[ 

Using this and the multiplicities of the eigenvalues on the Fermat quintic Q p given 
in ITable 2\ we can now perform the Z5 x Zs-quotient and obtain the degeneracies of 
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the scalar Laplacian on the Qp- The results are listed in ITable~6l We find complete 
agreement with the spectrum found by directly computing the eigenvalues on Qf 
given in ITable 4[ Naturally, this comparison is limited by the number of eigenvalues 
we were able to compute on Qp. The agreement of the lower lying levels, however, 
gives us confidence that the values of and /i^ for m = 3, 4, 5, 6, . . . given in ITable 4l 
are also a good approximation to the exact results on the quotient. 



5.5 Varying the Complex Structure 

To numerically compute any metric-dependent quantity on a Calabi-Yau manifold, one 
has to fix the complex structure and Kahler moduli to specific values. This was done, 
for example, in [Subsection 5.31 where the moduli were chosen so that the covering 
space was the Fermat quintic with unit volume. In this section, we will extend our 
results to the one-parameter family of Z5 x Z5 symmetric quintics defined by the 
vanishing of the polynomial 

The Kahler modulus will always be fixed so that the volume of Q^p is unity. The 
complex structure parameter if) can, in principle, take on any complex value. However, 
for simplicity, we will only consider -j/^ G M in this subsection. Note that each is 
indeed a quintic with the free Z5 x Z5 symmetry in eq. (176|) . Hence, the quotient 



Q^ = Q^I (Z5 X Z5) (105) 

is a smooth Calabi-Yau threefold. 

We have computed the spectrum of the scalar Laplace operator on this one- 
parameter family of quotients for various values of ^. The resulting ?/;-dependent 



spectrum can be found in Figure 13 Note that this one-parameter family of Z5 x Z5- 
symmetric quintics passes through two special points, 

'0 = 0: Without the term, (5^=o = is exactly the Fermat quintic. We will 
investigate the symmetry enhancement at this point in the next subsection. 

•0 = 1: This is the so-called conifold point, where the quintic is singular. On the 
covering space Q^=\ C P^, the singularity is at 

^C7 = [1 : 1 : 1 : 1 : 1] (106) 

and its images under the Z5 x Z5 symmetry group. At these points the hyper- 
surface equation fails to be transversal, 

^l^(^c) = ■ ■ ■ = ^l^(^c) = g^=i(^c) = 0, (107) 
causing the singularity. 
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Figure 13: Spectrum of the scalar Laplace operator on the real 1-parameter 
family of quintic quotients. The metric is computed at degree 
kh = 10 with Uh = 406,250. The Laplace operator is evaluated at 
kff, = 10 and = 50,000 points. 
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Perhaps surprisingly, the spectrum of the scalar Laplace operator shows no trace of 
the conifold singularity at ip = 1. However, the reason for this is straightforward. The 
low-lying modes are slowly-varying functions and, in particular, are almost constant 
near any point-like singularity. For example, the first massive eigenvalue is essentially 
determined by the diameter of the manifold, see [Subsection 7.2\ and does not depend 
on local details of the metric. 



5.6 Branching Rules 



Let us return to spectrum of the Laplace-Beltrami operator in Figure 13 and focus 
on the neighbourhood of ip = 0. Clearly, Q^=o = Qf is the quotient of the Fermat 
quintic, while Q^^o is a deformation of the Fermat quotient that breaks part of its 
discrete isometry group. In particular, note that for small non-zero values of ip 

• The first massive level ^i{Qf) = 4 splits into two pairs of eigenvalues. 

• The second massive level ^2{Qf) = 2 remains two-fold degenerate. 

In this subsection, we will attempt to understand this from the group-theoretical 
perspective. 

As discussed in [Subsection 5.4| the multiplicities of the eigenvalues on the quo- 
tient Q^, = (5vV(^5 X ^5) are really determined by the representation theory of the 
symmetry group of the covering space. We have to distinguish two cases. 

ip = 0: This is the case of the Fermat quintic, whose symmetries we already discussed 
in [Subsection 4.3 



Aut (Q^=o) = Aut {Qf) = {S5 x Z2) x {Z^Y- (108) 



The irreducible representations of Aut [Qf) were presented in [Table's! 

i/j ^ 0: In this case, the invariance of the Yl Zi monomial gives one further constraint 
on the (Zs)"^ phase rotations. In other words, turning on breaks the phase 
rotation symmetry to (^5)^. The remaining symmetry group i£^ 

AHt (Q^^o) = (^5 X Z2) X (Z5)'. (109) 

The irreducible representations of Aut (Q^^o) are given in [Table 71 Note that, 
by construction, this group is a proper subgroup of the full symmetry group, 
both of which containing the free Z5 x Z5 action. That is, 

Ant(g^,=o) 3 Ant(g^^o) 3 Z5XZ5. (110) 
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d 


1 


4 


5 


6 


20 


24 


30 


40 


48 


60 


rid 


4 


4 


4 


2 


8 


2 


8 


4 


2 


2 




1 





1 


2 





4 


2 








4 



Table 7: Number rid of distinct irreducible representations of Aut(Q^,^Q) in 
complex dimension d. We also list the dimension dim^'"'^^^ of the 
Z5 X Z^-invariant subspace for each representation. Note that it turns 
out to only depend on the dimension d of the representation. 



AutiQp) 3 Ant(Q^^o) Aut{QF) 3 Ant(g^^o) 

1 — ^ 1 20 — ^20 

2 — > 1©1 30 — ^30 

4 — > 4 40i — ^40 

5 — ^5 4O2 — ^20020 

6 — ^6 6O1 — ^60 

8 — ^4©4 6O2 — ^30030 

10 — > 5©5 80 — > 40040 

12 — > 6©6 120 — > 48©48©24 

Table 8: Branching rules for the decomposition of the irreducible representa- 
tions o/Aut((5i?) into the irreducible representations of its subgroup 
Aut(Q^^o)- Note that there are always numerous distinct representa- 
tions in each dimension, see \Table 51 and[^ In particular, in dimen- 
sion 60 there are 10 irreps of Aut{Qp), which we denote by 60 1, that 
remain irreducible under Ant{Q^^o) and 8 irreps, denoted by 6O2, 
that branch into two distinct 30 -dimensional irreducible representa- 
tions. 
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As one turns on the deformation, the eigenvalues must spht according to the group- 
theoretical branching rules. We list these in ITable 81 

Finally, we are really interested in the eigenvalues on the quotient Q^, which means 
that one must restrict to the Z5 x Zs-invariants of each representation. For the Fer- 
mat quintic, we listed the number and the dimension, dim^^^^'', of these invariants 
in ITable 51 We list the analogous information for the Z5 x Zs-invariants within the 
irreducible representations of Aut((5^^o) in ITable "71 This allows us to compute the 
splitting of the eigenvalues on the quotient Q^. However, just knowing the multiplici- 
ties turns out to be not quite enough since same- dimensional but different irreducible 
representations can branch in different ways. In particular, the first massive level 
on Q^=Q comes from a 60-dimensional representation of Q^=o, which can branch in 
two ways according to ITable 81 However, since we have seen in [Figure 13] that the 



eigenvalues do branch, this 60-dimensional representation must be of the type GOg. 

To summarize, these group theoretical considerations are completely compatible 
with the observed branching of the eigenvalues under the complex structure deforma- 
tion by The low-lying eigenvalues of the scalar Laplacian on split as 

Z5XZ5 ^ Ant(Q^=o) Ant(g^^o) 3 Z5XZ5 
invariants irreps irreps invariants 

Ato(Qi/>=o) = 1 C 1 ^1 D Ato(Qv7^o) = 1 



30 




/^i (Q^=o) = 4 C _ 

^30 D /ii(Qv'^o)=2 
/i2(Q^=o)=2 C 30 -30 D ^iz{Q^^^^) =2. 

5.7 Another Family 

Finally, let us consider another family of complex structure moduli. First, we deform 
the Fermat quintic to a generic Z5 x Z5 invariant polynomial; that is, switch on all 
coefficients in eq. f l76l) . Then restrict to the real one-parameter family of covering 
spaces defined by 

Qf= 5Z + n + w{zl^iZi + eye) 

+ (1 - i)^{zlzizl + eye) - (1 - 2i)^{zlzlz:i + eye) (112) 
- (2 - i)if[zlz2Z^ + eye) 

and form the quotient spaces 

= Q^l (^5 X Z5). (113) 



19 



Since we chose ^ to be real, the complex conjugation Z2 remains unbroken. 
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For generic values of v?, this breaks all symmetries of the Fermat quintic except for 
the free Z5 x Z5 that we are dividing out. Consequently, we expect no degeneracies in 
the spectrum of the Laplace-Beltrami operator. In Figure 14, we plot the spectrum of 
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Figure I4: Spectrum of the scalar Laplace operator on the real 1-parameter 
family of quintic quotients. The metric is computed at degree 
kh = 10, n/i = 406,250 and the Laplace operator evaluated at = 
10 and = 50,000 points. 

A and, indeed, observe that the degeneracies of the eigenvalues on the Fermat quintic 
Qip=Q are broken as ip is turned on. 

6 A Heterotic Standard Model Manifold 

In this last section, we will compute the spectrum of the Laplace-Beltrami operator on 
the torus- fibered Calabi-Yau threefold X with vri(X) ='L^^'L^ that was used in |50] 
to construct a heterotic standard model. The threefold X is most easily described in 
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terms of its universal cover X, which is the complete intersection 

X = {P(x,t,y) = = C X Pf,^,^] X ^l,.,y,.,y,^ (114) 

defined by the degree-(3, 1, 0) and (0, 1, 3) polynomials 

P{x, t, y) = to + xl + xl + A1X0X1X2 j + tiAs (^xlx2 + xfxo + xlxi 

R{x, t, y) = ti (^yQ + yi + y2 + ^2ymy2) + to (vlvi + yly2 + ylyo) ■ 



(115) 



Note that Ai, A2, A3 G C end up parametrizing the complex structure of X. For 
generic Aj, the two maps 

{[xo : Xi : X2] ^ [xo : ooxi : uj'^X2] 
[to : ti] ^ [to : ujti] (116a) 
[vo ■ yi ■ 2/2] ^ [yo ■■ ^yi ■ ^^^^2] 

and 

{[xo : Xi : X2] ^ [xi : X2 : Xq] 
[to-.h] ^ [to-.h] (116b) 
[yo ■■ vi ■■ y2] ^ [yi ■ y2 ■ yo] 

generate a free Z3 x Z3 group action on X. Hence, the quotient 

X = X/(Z3XZ3) (117) 

is a smooth Calabi-Yau threefold. In addition to the /i^'^(X) = 3 complex structure 
moduli of X, there are also h^'^{X) = 3 Kahler moduli. The Kahler class on the alge- 
braic variety X is determined by a line bundle £ whose first Chern class is represented 
by the Kahler class, 

ci(£) = [ujx] G H^'\X, Z) = H^'\X, C) n H^{X, Z). (118) 

Pulling back to the covering space with the quotient map q, the Kahler class is equiv- 
alently encoded by an equivariantl^ line bundle 

q*{L) =0^(ai,6,a2), (119) 

which is determined by some ai, b, 02 G Z>o- Note that, by definition, the sections of 
0^(ai, b, 02) are the homogeneous polynomials in x, t, and y of multidegree (ai, b, 02)- 
We now want to compute the Calabi-Yau metric on the quotient X using Don- 
aldson's algorithm. However, as discussed in detail in the previous section, we will 



^''As was shown in [42l 140] , equivariance requires ai + 02 = mod 3. We will always use the 
equivariant action specified by eqns. pi6ap and (jll6b|) . 
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formulate everything in terms of Z3 x Zs-invariant data on the covering space X. 
First, one has to pick a muhidegree 

kh = {ci'i,b,a2) G (Z>o)'^7 ai + 02 = mod 3 (120) 
determining the Kahler class of the metric. Then one has to find a basis 

span {sQ,|a = 0, . . . , N{kh) — l} = 

= C[xo, xi, X2, to, h, yo, yi, l/2]fcf {R{x, t, y), P{x, t, y)) (121) 

for the invariant sections of 0^(ai, 6, 02) modulo the complete intersection equations, 
as described in detail in [IQ]. This is all the data needed to apply Donaldson's algo- 
rithm and compute the approximate Calabi-Yau metric. Note that, since we always 
normalize the volume to unity, the exact Calabi-Yau metric only depends on the ray 
Qkh but not on the "radial" distance gcd(ai, 6, 02). However, the number of sections 
N{kh) and, therefore, the number of parameters in the matrix does depend on 
kh explicitly. Going from kh to 2kh, 3kh, ■ ■ ■ increases the number of parameters 
and subsequently improves the accuracy of the Calabi-Yau metric computed through 
Donaldson's algorithm. 



6.1 The Spectrum of the Laplacian on X 

Having determined the metric, we now turn towards the spectrum of the Laplace- 
Beltrami operator. We do this again by computing the matrix elements of the 
Laplacian on the covering in an approximate basis of Z3 x Zs-invariant functions, 
completely analogous to [Subsection 5.11 To specify the truncated space of invariant 
functions on X, fix a multidegree k^ proportional to kh] that is, 

k^ = (k^,, A;^2' ^03) e '^kh n (Z>o)'. (122) 

Then pick a basis {saja = 0, . . . ,N{kfj,) — l} of degree-Zc,^ homogeneous, Z3 x Z3- 
invariant polynomials. These define a finite-dimensional space of invariant functions 
on X as 



Xi 



a,P = 0,...,N{k^)-l} . (123) 



By computing the matrix elements of the Laplacian and solving the (generalized) 
matrix eigenvalue problem, we obtain the eigenvalues A^^^^^* of the Laplacian on the 
covering space X acting on Z3 x Z3-invariant functions. These are identical to the 
eigenvalues of the Laplacian on X, but with volume 

^°"^)=E7^''°'<-^>- ''''' 
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In the computation on X we normalized the volume to unity. Hence, after rescaling 
the volume of X back to one, the eigenvalues of the scalar Laplacian on X are 

\Z3xZ3 

A,.^^. (125) 



In Figure 15 , we compute the spectrum of the Laplace-Beltrami operator on X at 
two different points in the Kahler moduli space but with the same complex structure. 
Recall that we always normalize the volume, corresponding to the "radial" distance in 
the Kahler moduli space, to unity. The non-trivial Kahler moduli are the "angular" 
directions in the Kahler cone, and we consider the two different rays Q ■ (2, 1, 1) and 
Q ■ (2, 2, 1). As expected, the actual eigenvalues do depend on the Kahler moduli, as 



is evident from Figure 15 



Furthermore, note that there appear to be no multiplicities in the spectrum. At 
first sight, this might be a surprise to the cognoscente, as there is a residual symme- 
try. By construction [l2], the covering space X comes with a (^3)^ group action of 
which only a Z3 x Z3 subgroup acts freely and can be divided out to obtain X. The 
remaining generators are 

{[xq : Xi : X2] ^ [xi : X2 : Xq] 
[to:ti] ^ [to-ti] (126a) 
[vo ■■ vi ■ y2] ^ [yo ■■ yi ■ y2] 

and 

{[xo : Xi : X2] ^ [xo : Xi : X2] 
[to-.ti] ^ [to-.ti] (126b) 
[yo ■■ yi ■ 1/2] ^ [yi ■ 1/2 : yo] 

in addition to 71 and 72, see eqns. flll6ap and (1116b[) . Moreover, we used the point 
Ai = = A2, A3 = 1 where the polynomials eq. flll5p are also invariant under complex 
conjugation. Hence, the symmetry group on the covering space is 

'MtiX) = Z2 i< {ZsY = Dg X {ZsY . (127) 

To understand the latter identity, note the Z2 action in the semidirect product: 

• Complex conjugation commutes with 72, 73, and 74. 

• Complex conjugation does not commute with 71, but satisfies 



71 (^[xo : xi : X2], [to : h], [yo, ^2] j = ll[[xo ■ Xi : X2], [to : ti], bo, ^1,2/2] 

(128) 

Hence, 71 together with complex conjugation generate Dq, the dihedral group 
with 6 elements. 
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Figure 15: Eigenvalues of the scalar Laplace operator on the Z3 x Z^-threefold 
X with complex structure Ai = = A2, A3 = 1 and at two distinct 
points in the Kdhler moduli space. The metric is computed at degree 
kfi = (6,3,3) and = 170,560 points as well as degree = 
(6, 6, 3) and = 290,440 points, corresponding to the two different 
Kdhler moduli. The matrix elements of the scalar Laplacian are 
always evaluated on = 500,000 points. The blue pluses and 
crosses, corresponding in each case to k^ with the largest radial 
distance, are the highest precision eigenvalues for the two metrics. 



Aut(X)-Rep. 



Pl,---,P36 P37,---,P54 P55, 



.P81 



dim(p) 
dim (p^3xZ3^ 



1 




2 




Table 9: Number of distinct irreducible representations of Aui{X) in com- 
plex dimension d. We also list the dimension dim^^^^^ of the Z3 XZ3- 
invariant suhspace for each representation. 
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The group Aut(X) is of order 162 = 6 x 3'^ and has one- and two-dimensional repre- 
sentations due to the Dq factor. As discussed previously, the surviving eigenfunctions 
on the quotient X are the Z3 x Zs-invariant eigenfunctions on the covering space X. 
Hence, we have to determine the subspace invariant under the freely acting Z3 x Z3 
inside of Aut(X). We list all this data in ITable"9l We find that all the multiplicities 
on X are, indeed, one. 



7 The Sound of Space-Time 

7.1 Kaluza-Klein Modes of the Graviton 

Consider a 10-dimensional spacetime of the form M^'^ x Y, where Y is some real, 
compact 6-dimensional Calabi-Yau manifold. Since Y is compact, there is a scale 
associated with it. Let us agree on a unit of length L such that Vol(y) = 1 ■ L^. The 
gravitational interactions in this world are complicated, but have two easy limiting 
cases. First, if the separation r of two probe masses Mi and M2 is large, then the 
gravitational potential between them is given by Newton's law 

» i) . _G,Mi^§. (129) 

r 

In the other extreme, when r is very small, the potential becomes the Green- Schwarz- 
Witten law 

V{r^L) = -Gi,^^. (130) 

By dimensional analysis 

G4 ~ ^ , (131) 

with a constant of proportionality independent of Y to be determined below. In- 
between these two extremal limits for the separation r, the gravitational potential is 
a complicated interpolation between eq. fll29p and eq. (11301) . 

There are two alternative ways of describing fields on M^'^ x Y. One can either 
directly use 10-dimensional field theory, or work with an infinite tower of massive 
Kaluza-Klein fields depending on M'^'^ only. Both methods are equivalent, but for the 
purposes of this paper we only consider the Kaluza-Klein compactification [STl [521153] . 



In this approach, the single 10-dimensional massless graviton g^^^\ A, 5 = 0, . . . , 9 
is decomposed into 4-dimensional gravitons, vectors, and scalars. For simplicity, let 
us only consider 4-dimensional gravity, that is, 4-d fields with symmetrized indices 
a, 6 = 0, . . . , 3. Then 

00 

gT' (-0, . . . , X3, yi, . . . , ye) = E • • • ' ^6) ■ ^,"^'"(xo, . . . , xs), (132) 

n=0 
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where the . . . , ye) G 1^-dependence of the 10-dimensional metric is now encoded 
in a basis of functions G C°°(y, M). The most useful such basis consists of the 
solutions to the equations of motion on Y , that is, the eigenf unctions of the scalar 
Laplace operator 

Ay0„(?/i, . . . ,?/6) = A„0„(?/i, . . . ,?/6), A„ < A„+i. (133) 
The corresponding 4-dimensional Lagrangian contains the infinite tower of fields 

(4D),n f 

g^^ij or mass 

ITT'n = n = 0, . . . , oo. (134) 

As discussed previously, there is a unique zero mode Aq = leading to a single massless 
graviton in 4 dimensions. The gravitational potential is then the sum of the potential 
due to the massless graviton plus the Yukawa-interaction of the massive modes. 



V(r) = -G,^±e-"^' = -g/Ji^ i 1 + te-'"'A ^ (135) 

n=0 \ n=l / 

At distance scales ^ only the massless graviton propagates. This expected 
behaviour is clearly visible in the r ^ limit of eq. fll35p . and one immediately 
recovers eq. f ll29p . At distance scales r <^ on the other hand, the massless graviton 
as well as the infinite tower of massive spin-2 fields propagate. The corresponding 
asymptotic behaviour of the gravitational potential is less obvious. However, note 
that the asymptotic growth 

lim = 384:71^1-^ ^ m.,, > 2^fl^n^l^L-^ (136) 

of the Kaluza-Klein masses is known from Weyl's formula, see [Subsection 3."3l Hence, 
the r <^ ^ limit of eq. fll35p is 

n=0 



^ r Jn=o ^ Stt^ 




Again, this matches the expected behaviour eq. (11301) . 

The purpose of this section is to fill the gap between the extremal limits and 
determine the gravitational potential at distances r ^ L. This explicitly depends 
on the details of the internal Calabi-Yau threefold Y, and there is no way around 
solving eq. (I133p . The eigenvalues A„ and corresponding eigenfunctions 0„ depend 
on the Calabi-Yau metric and can only be computed numerically. We have presented 
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Figure 16: The gravitational potential V {r) computed from eq. (|135|) on. 



d3,1 



Qf, where Qp is the Fermat quintic with unit volume, YoI{Qf) = 
1 • L^. The Kaluza-Klein masses run = y/X^ are computed using 



the numerical results for A„ given in Subsection 4-2 
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a detailed algorithm for calculating the spectrum of A in this paper, and given the 
results for a number of different Calabi-Yau threefolds. As an example, let us compute 
the gravitational potential V{r) derived from the numerical eigenvalues of the scalar 
Laplace operator on the Fermat quintic discussed in [Subsection 4.2[ The result is 



plotted in Figure 16 



7.2 Spectral Gap 



As is evident from [Figure 16 deviations from the pure ^ (green line) and ^ (red line) 



potentials occur for r in the region where these gravitational potentials have a similar 
magnitude. In fact, these curves intersect at 



^ M1M2 IhGiL'' M1M2 6/15 , , 

Note that this point of intersection is independent of the Calabi-Yau manifold and its 
geometry. As will become clear below, for Calabi-Yau threefolds which are relatively 
"round" , such as the Fermat quintic, tq is a good estimate for the point of substantial 
deviation from the ^ potential. However, for geometries that are stretched or develop 
a throat in at least one direction, this deviation point is best determined by another 
scale, in principle independent of the volume of the internal space. This other scale is 
the mass mi of the lightest Kaluza-Klein modj^. see eq. fll35p . For such manifolds, 
the spectral gapE§ Ai and, hence, the mass mi becomes smaller. Eventually, the 
manifold may be sufficiently elongated that 3> rg. In this case ^ becomes the 
best estimate of the point of deviation from the ^ potential. 

Of course, both the volume and Ai = m\ are determined by the geometry of the 
internal Calabi-Yau manifold. However, what geometric property really determines 



^^Thc leading order correction to the gravitational potential is often [5^ parametrized by the 
lowest Kaluza-Klein mass mi and its multiplicity as 

y(r)«-G4^^^^(l + Mie"'"'0- (139) 

While this works well for symmetric spaces like spheres and tori with their large multiplicities and 
widely-separated eigenvalues, there are two issues when dealing with more general manifolds: 

• The multiplicity is caused by symmetries, and tiny non-symmetric deformations can (and will) 
make the eigenvalues non-degenerate (sec lScction b\ . 

• The separation between the zero mode and the first massive mode is, in general, much larger 
than the separation between the first and second mode. For example, on the non-symmetric 
"random quintic" Calabi-Yau threefold in [Subsection 4.H 

7710 = 0, mi « 5.95, ma « 6.00. (f40) 

^^The first massive eigenvalue of the scalar Laplacian, Ai, is also called the spectral gap since it is 
the gap between the unique zero mode Aq = and the first massive mode. 
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the spectral gap Ai? In fact, this is determined by the "diameter" of the mani- 
fold. Recall that the diameter D is defined to be the largest separation of any two 
points, as measured by the shortest geodesic between them. Then, on an arbitrary 
real ci- dimensional manifold with non-negative scalar curvaturj^. the spectral gap is 
essentially determined by the diameter via [55l [56| [57] 

vr^ ^ 2d{d + A) tx J2d(d + 4) , , , 

Clearly, in a compactification where all internal directions are essentially of equal size, 
the diameter is of the order of 1 ■ L. However, as soon as there is even one elongated 
internal direction or one long throat/spike develops, the diameter can be very large. 
Hence, the spectral gap becomes very small and deviations from ^ gravity appear for 
relatively large values of r ~ 

The definition of the diameter D is very impractical if one wants to explicitly 
calculate it, since this would require global knowledge about the shortest geodesies. 
However, to get a rough estimate of D, one can reverse the inequalities eq. fll4ip 
and then use the numerically computed value for Ai. For example, on the Fermat 
quintic our numerical computation in [Subsection 4.3] yielded Ai ~ 41.1. Therefore, 
the diameter must be in the range 

0.490 ^^<D< ^^'-'^^'^ ^ 1.71. (142) 
vAi vAi 

Thus, computing the value of Ai numerically on a Calabi-Yau threefold for specific 
values of its moduli gives us direct information about the "shape" of the manifold; 
information that would be hard to obtain by direct calculation of the diameter D. 
For example, it follows from eq. (11421) that the Fermat quintic is relatively "round" . 
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A Spectrum of the Laplacian on Projective Space 

In this Appendix, we compute the lowest eigenvalue of the Laplace operator on 
using the rescaled Fubini-Study Kahler potential eq. f[T6l) . To do this, go to the 

^^In particular, a Calabi-Yau |-fold. 
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coordinate patch where zq = 1 and use zi, Z2, z^ as local coordinates. We find that 

'l + l^ip 2:2^1 2:3^1 

^^i = ^7r(l + |zi|2+|^2|' + |z3n I ^21^2 1 + ^21' Z3Z2 



zizs Z2Z3 1 + \z3\''J (143) 
det (oi^) = -. — 

and, hence, 

A = (d-,g-'^ det{g)d, + d.g'^^ det{g)d^~) . (144) 

One can now compute the eigenvalue corresponding to the eigenfunction 0i i in 
eq. ([27D. We find that 



A0i,i = 2—— (d,g'^ det{g)di + d,g^^ det{g)d-^ - 



|2 



-z, ^^^^^ 



1 + |2;i|2 + l^sP + |2;3|2' 
Hence, is indeed an eigenfunction of A with eigenvalue 

IGtt 47r , , , , 

Hence, the numerical coefficient in eq. (|T8i) is indeed the correct one for our volume 
normalization Voli^(P^) = 1. 



B Semidirect Products 

Let G and be two groups, and let 

tlj-.G^ Aut(iV) (147) 
be a map from G to the automorphisms of A^. The semi-direct product 

G N = \^{n,g) neN,geG^ (148) 

is defined to be the group consisting of pairs (ra, g) with the group action 

{ni.gi) ■ {n2,g2) = {rii ■ ^Pigi){n2), gi ■ 5-2)- (149) 

Usually, one just writes G N with the map implied but not explicitly named. 
Note that G is a subgroup and is a normal subgroup of the semidirect product. 
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For example, consider the semidirect product with = 3^ and = (^5)^ used in 
[Subsection 4.3[ These two groups are acting on five homogeneous via permutation^ 
and phase rotations 

/ \ i- 2Tvin-^ 2'Kin'2 27vin^ 27rin^ _. 

(^(ni,n2, rig, 71,4), [zq, Zi, Z2, z^, z^]j ^ [zq, Zie^ , Z2e^ , z^e^ , Z4e^ \, (150) 

respectively. The two group actions do not commute, and, therefore, the total sym- 
metry group is not simply the product S5 x (Zs)^. The "non-commutativity" between 
5*5 and (1^5)^ is encoded in a map 

ip : ^ Aut ^(Zs)^^ , a (^n a^^ o n o . (151) 

To be completely explicit, note that the permutation group is generated by the 
cyclic permutation c and a transposition t, acting as 

t- [zq, Zi, Z2, Z3, Zi'\ I— i> [zq, Zi, Z2, Z4^, z^j, (152) 
C: [2:0, ^1, ^2; ^3; ^4] ^ [Zl, Z2, Zs, Zi, Zq^. 

The generators {c,t) = S5 act, via ip, on (Zs)"^ as 

ipit) : (Zg)^ (Zg)^ (ni, n2, rig, n^) (ni, n2, n4, n^) ^^^^^ 
ip{c) : {Z^y (Zs)^, (ni, n2, ris, n4) t-^ (-^4, ni - n^, n2 - ^4, ns - n4) 

It is straightforward, if tedious, to show that is a. group homomorphism and that 
the total symmetry group generated by 5*5 and (Zs)^ is, in fact, the semidirect product 

55 {Z,)\ (154) 

By the usual abuse of notation, we always drop the subscript ip in the main part of 
this paper. 

C Notes on Donaldson's Algorithm on Quotients 

For explicitness, let us consider the same setup as in [Subsection 5.H that is, Q C 
is a Z5 X Z5 symmetric quintic and we want to compute the metric on the quotient 
Q = Q/ (Z5 X Z5). To fix notation, let us denote the two generators for the character 
ring of the group by 

Xl(^7l)=e^'^^/^ Xi{92) = l, 

X2{gi) = l, X2{g2) = e^-''\ 

We consider homogeneous polynomials in degrees kh G 5Z, so there is a linear Z5 x Z5 
group action. In eq. 0821) we determined the invariant polynomials. Now, let us 



24 



S'5 is, by definition, the group of permutations of five objects. 
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slightly generalize this result and determine "covariant polynomials" transforming as 
some character x of the group, 

P ° 9{z) = X{9)p{z) geZ^xZ^. (156) 

These again form a linear space of ^-covariant polynomials, which we denote as 

C[^o, zi, Z2, Z3, Zi]\ = p o g{z) = xi9)p{z)y (157) 

Note that the covariant polynomials do not form a ring, but rather a module over the 
invariant ring. Nevertheless, by a slight generalization of the Hironaka decomposition, 
we can express the covariants clS 8b direct sum 

100 

C[zo, Zi, Z2, 23, Z^ll^ = C[6'i, 6*2, 6*3, 6*4, 6'5]fc^_deg(;,X), (158) 

1=1 

where the 9i, . . . , 9^ G C[zo, zi, Z2, z^, 24]^^^^^ can be taken to be the primary invari- 
ants of the original Hironaka decomposition eq. fl82l) and the "secondary covariants" 
7]^, . . . , ?7^QQ are certain x-covariant polynomials that need to be computed [58] • For 
example, we find 



Xl 4 I 4 I 4 I 4 I 4 

Xl 3 3 3 3 3 

7/2 = ZqZ-^Zs + Z1Z2Z4 + 2:2^3^0 + ^3^4^1 + ^4^o^2; • • • 



(159) 



and 



^71 = ^0 + e 5 + 5 ^2 + 5 23 + 5 2;^, 
?7^^ = 20-21-22 + e~ 21^2^3 + e "2:22:324 + e "232420 + e ~ 242521, 



(160) 



Note that we always take the defining quintic polynomial Q{z) to be completeljo 
invariant, see eq. (!76l) . Restricting everything to the hypersurface Q{z) = 0, we get 
homogeneous polynomials on the Calabi-Yau threefold. We pick bases {s^} for the 
X-covariant polynomials, that is, 

X = 1 : span {si} = C[zo, Zi, 22, 23, 24]^^""^^ {Q{z)), 

XT^l: span {s^} = (^C[2o, 21, 22,23, 24]fc;,/(Q(2)))^ (^61) 

= C[2o, 2i, 22, 23, 24]^^. 



^^If Q{z) were a ^-covariant polynomial, it would still define a Z5 x Z5 invariant Calabi-Yau hy- 
persurface. Everything in this paper would generalize straightforwardly, so we ignore this possibility 
to simplify notation. 
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We now turn towards computing the metric on the quotient Q or, equivalently, 
computing the Z5 x Zs-invariant metric on the covering space Q by a variant of 
Donaldson's algorithm. For this, we pick the ansatz 

4 4 

. xixt 

^(^'^) = - E (162) 

X=xU°2 

for the Calabi-Yau metric. One can think of /i as a block- diagonal matrix with blocks 
labelled by the characters x- The T-operator is likewise block-diagonal, and therefore 
one obtains a balanced metric as the fixed point of the iteration 

hr'^hrJ = T{hrT'- (163) 

Note that this fixed point is the samj^ as what one would obtain from Donaldson's 
algorithm on the covering space Q (without using any symmetry). Only now the basis 
of sections is such that the impact of the Z5 x Z5 symmetry is clearly visible: h is 
block- diagonal with blocks labelled by the characters x- 

As usual, the balanced metrics are better and better approximations to the Calabi- 
Yau metric as one increases the degree k^. We find that this method of computing 
the Calabi-Yau metric on the quotient Q is the most effective. 



Bibliography 

[1] p. Candelas, G. T. Horowitz, A. Strominger, and E. Witten, "Vacuum 
Configurations for Superstrings," Nucl. Phys. B258 (1985) 46-74. 

[2] A. Lukas, B. A. Ovrut, and D. Waldram, "On the four-dimensional effective 
action of strongly coupled heterotic string theory," Nucl. Phys. B532 (1998) 
43-82, hep-th/9710208 , 

[3] R. Donagi, A. Lukas, B. A. Ovrut, and D. Waldram, "Non-perturbative vacua 
and particle physics in M-theory," JHEP 05 (1999) 018, hep-th/9811168t 

[4] A. Lukas, B. A. Ovrut, and D. Waldram, "Non-standard embedding and 
five-branes in heterotic M- theory," Phys. Rev. D59 (1999) 106005, 



ihep-th/9808101 



[5] A. Lukas, B. A. Ovrut, K. S. Stelle, and D. Waldram, "The universe as a 



domain wall," Phys. Rev. D59 (1999) 086001, hep-th/9803235 



^^And different from the fixed point where one restricts to only the invariant sections. The latter 
is just the X = 1 block. 



62 



[6] A. Lukas, B. A. Ovrut, and D. Waldram, "The ten-dimensional effective action 
of strongly coupled heterotic string theory," Nucl. Phys. B540 (1999) 230-246, 



hep-th/9801087 



[7] R. Donagi, A. Lukas, B. A. Ovrut, and D. Waldram, "Holomorphic vector 
bundles and non-perturbative vacua in M- theory," JHEP 06 (1999) 034, 



hep-th/9901009 



R. Donagi, B. A. Ovrut, and D. Waldram, "Moduli spaces of fivebranes on 



elliptic Calabi-Yau threefolds," JHEP 11 (1999) 030, hep-th/9904054 



[9] R. Donagi, B. A. Ovrut, T. Pantev, and D. Waldram, "Standard models from 
heterotic M-theor y," Adv. Theor. Math. Phys. 5 (2002) 93-137, 
hep-th/9912208[ 



[10] R. Donagi, B. A. Ovrut, T. Pantev, and D. Waldram, "Standard-model 
bundles," Adv. Theor. Math. Phys. 5 (2002) 563-615, linath/OOOBOlOl 



[11] E. I. Buchbinder, R. Donagi, and B. A. Ovrut, "Vector bundle moduli 

superpotentials in heterotic superstrings and M-theory," JHEP 07 (2002) 066, 
]hep-th/0206203| 

[12] R. Donagi, B. A. Ovrut, T. Pantev, and D. Waldram, "Standard-model bundles 
on non-simply connected Calabi-Yau threefolds," JHEP 08 (2001) 053, 
hep-th/0008 008[ 



[13] R. Donagi, B. A. Ovrut, T. Pantev, and D. Waldram, "Spectral involutions on 
rational elliptic surfaces," Adv. Theor. Math. Phys. 5 (2002) 499-561, 



math/00080Tll 



[14] B. A. Ovrut, T. Pantev, and R. Reinbacher, "Invariant homology on standard 



model manifolds," JHEP 01 (2004) 059, hep-th/0303020 



[15] R. Donagi, Y.-H. He, B. A. Ovrut, and R. Reinbacher, "The particle spectrum 
of heterotic compactifications," JHEP 12 (2004) 054, hep-th/0405014 . 

[16] R. Donagi, Y.-H. He, B. A. Ovrut, and R. Reinbacher, "Moduli dependent 
spectra of heterotic compactifications," Phys. Lett. B598 (2004) 279-284, 



hep-th/0403291, 



[17] R. Donagi, Y.-H. He, B. A. Ovrut, and R. Reinbacher, "The spectra of 



heterotic standard model vacua," JHEP 06 (2005) 070, .hep-th/0411156 



[18] R. Donagi, Y.-H. He, B. A. Ovrut, and R. Reinbacher, "Higgs doublets, split 
muhiplets and heterotic SU(3)C x SU(2)L x U(1)Y spectra," Phys. Lett. B618 



(2005) 259-264, hep-th/0409291, 



63 



V. Braun, Y.-H. He, B. A. Ovrut, and T. Pantev, "Heterotic standard model 
moduli," JHEP 01 (2006) 025, hep-th/0509051 . 

V. Bouchard and R. Donagi, "An SU(5) heterotic standard model," Phys. Lett. 
B633 (2006) 783-791, .hep-th/0512149, 



V. Braun, Y.-H. He, B. A. Ovrut, and T. Pantev, "A standard model from the 
E(8) X E(8) heterotic superstring," JHEP 06 (2005) 039, ,hep-th/0502155| 

V. Braun, Y.-H. He, B. A. Ovrut, and T. Pantev, "The exact MSSM spectrum 



from string theory," JHEP 05 (2006) 043, |hep-th/0512177 



V. Braun, Y.-H. He, B. A. Ovrut, and T. Pantev, "A heterotic standard 



model," Phys. Lett. B618 (2005) 252-258, [hep-th/0501070 



P. Candelas and S. Kalara, "Yukawa couplings for a three generation 
superstring compactification," Nucl. Phys. B298 (1988) 357. 

V. Braun, Y.-H. He, and B. A. Ovrut, "Yukawa couplings in heterotic standard 



models," JHEP 04 (2006) 019, |hep-th/060T204 



V. Braun, Y.-H. He, B. A. Ovrut, and T. Pantev, "Moduli dependent mu-terms 
in a heterotic standard model," JHEP 03 (2006) 006, ,hep-th/0510142, 

V. Bouchard, M. Cvetic, and R. Donagi, "Tri-linear couplings in an heterotic 
minimal supersymmetric standard model," ,hep-th/ 0602096^ 

P. Candelas, X. C. De La Ossa, P. S. Green, and L. Parkes, "A pair of 
Calabi-Yau manifolds as an exactly soluble superconformal theory," Nucl. Phys. 
B359 (1991) 21-74. 

B. R. Greene, D. R. Morrison, and M. R. Plesser, "Mirror manifolds in higher 
dimension," Commun. Math. Phys. 173 (1995) 559-598, hep-th/9402ll9[ 

R. Donagi, R. Reinbacher, and S.-T. Yau, "Yukawa couplings on quintic 



threefolds," hep-th/0605203^, 



S. K. Donaldson, "Some numerical results in complex differential geometry," 
|iath.DG/05 126251 

M. R. Douglas, R. L. Karp, S. Lukic, and R. Reinbacher, "Numerical solution 



to the hermitian Yang-Mills equation on the Fermat quintic," hep-th/0606261 



M. R. Douglas, R. L. Karp, S. Lukic, and R. Reinbacher, "Numerical 
Calabi-Yau metrics," ,hep-th/0612075, 



64 



[34] S. K. Donaldson, "Scalar curvature and projective embeddings. II," Q. J. Math. 
56 (2005), no. 3, 345-356. 

[35] G. Tian, "On a set of polarized Kahler metrics on algebraic manifolds," J. 
Differential Geom. 32 (1990), no. 1, 99-130. 

[36] M. Headrick and T. Wiseman, "Numerical Ricci-flat metrics on K3," Class. 
Quant. Grav. 22 (2005) 4931-4960, |hep-th/0506129( 

[37] C. Doran, M. Headrick, C. P. Herzog, J. Kantor, and T. Wiseman, "Numerical 
Kaehler-Einstein metric on the third del Pezzo," hep -th/ 0703057[ 

[38] X. Wang, "Canonical metrics on stable vector bundles," Comm. Anal. Geom. 
13 (2005), no. 2, 253-285. 

[39] B. Sturmfels, Algorithms in invariant theory. Texts and Monographs in 
Symbolic Computation. Springer- Verlag, Vienna, 1993. 

[40] V. Braun, T. Brelidze, M. R. Douglas, and B. A. Ovrut, "Calabi-Yau Metrics 



for Quotients and Complete Intersections," arXiv: 0712. 3563 [hep-th] 



[41] M. B. Green, J. H. Schwarz, and E. Witten, "Superstring Theory. Vol. 2: Loop 
Amplitudes, Anomalies and Phenomenology,". Cambridge, Uk: Univ. Pr. 
(1987) 596 P. (Cambridge Monographs On Mathematical Physics). 

[42] V. Braun, B. A. Ovrut, T. Pantev, and R. Reinbacher, "Elliptic Calabi-Yau 
threefolds with Z(3) x Z(3) Wilson lines," JHEP 12 (2004) 062, 
|he p-th/0410055[ 



[43] P. Candelas, X. de la Ossa, Y.-H. He, and B. Szendroi, "Triadophilia: A Special 



Corner in the Landscape," arXiv: 0706 . 3134 [hep-th] 



[44] C. luliu-Lazaroiu, D. McNamee, and C. Saemann, "Generalized Berezin 
quantization, Bergman metrics and fuzzy Laplacians," ,0804.4555,. 

[45] A. Ikeda and Y. Taniguchi, "Spectra and eigenforms of the Laplacian on 5*" and 
P"(C)," Osaka J. Math. 15 (1978), no. 3, 515-546. 

[46] E. Gabriel, G. E. Fagg, G. Bosilca, T. Angskun, J. J. Dongarra, J. M. Squyres, 
V. Sahay, P. Kambadur, B. Barrett, A. Lumsdaine, R. H. Castain, D. J. Daniel, 
R. L. Graham, and T. S. Woodall, "Open MPI: Goals, Concept, and Design of 
a Next Generation MPI Implementation," in Proceedings, 11th European 
PVM/MPI Users' Group Meeting, pp. 97-104. Budapest, Hungary, September, 
2004. 



65 



[47] E. Anderson, Z. Bai, C. Bischof, S. Blackford, J. Demmel, J. Dongarra, 

J. Du Croz, A. Greenbaum, S. Hammarling, A. McKenney, and D. Sorensen, 
LAPACK Users' Guide. Society for Industrial and Applied Mathematics, 
Philadelphia, PA, third ed., 1999. 

S. K. Donaldson, "Scalar curvature and projective embeddings. I," J. 
Differential Geom. 59 (2001), no. 3, 479-522. 

A. Kehagias and K. Sfetsos, "Deviations from the 1/r^ Newton law due to extra 
dimensions," Phys. Lett. B472 (2000) 39-44, lhep-ph/9905417| 

B. A. Ovrut, "A heterotic standard model," AIP Gonf. Proc. 805 (2006) 
236-239. 

T. Kaluza, "On the Problem of Unity in Physics," Sitzungsher. Preuss. Akad. 
Wiss. Berlm (Math. Phys. ) 1921 (1921) 966-972. 

O. Klein, "Quantum theory and five-dimensional theory of relativity," Z. Phys. 
37 (1926) 895-906. 

N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, "Phenomenology, 
astrophysics and cosmology of theories with sub-millimeter dimensions and TeV 
scale quantum gravity," Phys. Rev. D59 (1999) 086004, |hep-ph/9807344| 

F. Leblond, "Geometry of large extra dimensions versus graviton emission," 



Phys. Rev. D64 (2001) 045016, hep-ph/0104273 



J. Q. Zhong and H. C. Yang, "On the estimate of the first eigenvalue of a 
compact Riemannian manifold," Sci. Sinica Ser. A 27 (1984), no. 12, 
1265-1273. 

S. Y. Cheng, "Eigenvalue comparison theorems and its geometric applications," 
Math. Z. 143 (1975), no. 3, 289-297. 

M. Berger, A panoramic view of Riemannian geometry. Springer- Verlag, Berlin, 
2003. 

G.-M. Greuel, V. Levandovskyy, and H. Schonemann, "Singular: :Plural 
2.1," A Computer Algebra System for Noncommutative Polynomial Algebras, 
Centre for Computer Algebra, University of Kaiserslautern, 2003. 
http : //www. singular .uni-kl . de/plural. 



66 



